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SUMMARY 


This  report  summarizes  the  hover  performance  studies  for  the 
Advanced  Tandem  Helicopters  ranging  in  gross  weight  from  Medium 
Lift  Helicopters  (MLH)  of  50,000  pounds  gross  weight  to  Very 
High  Lift  Helicopters  (VHLH)  at  200,000  pounds,  and  describes 
the  full  flight  envelope  simulator  math  model  for  piloted  sim¬ 
ulation  of  these  configurations.  The  programming  and  the  model 
validation  approaches  are  also  discussed. 

Four  crane-type  tandem  rotor  helicopters  with  the  gross  weight 
capability  of  50,000  pounds,  80,000  pounds,  120,000  pounds, 
and  200,000  pounds  are  defined  in  terms  of  rotor  and  fuselage 
parameters  and  the  resulting  dynamic  characteristics.  An  auto¬ 
matic  Precision  Hover  System  (PHS)  is  synthesized  and  the  hover 
hold  system  performance  established  analytically  including 
aircraft  position  errors  and  pilot  station  acceleration  envir¬ 
onment  as  a  function  of  low  speed  turbulence.  In  addition, 
pilot  modeling  techniques  are  utilized  to  analytically  estimate 
the  piloted  performance  and  rating  for  various  attitude  and 
linear  velocity  control  response  modes  without  the  aid  of 
automatic  hold  functions.  These  estimated  performance  and 
rating  levels  are  compared  with  similar  data  obtained  from  a 
piloted  simulator  evaluation  on  the  Northrop  Corporation  (LAS/ 
WAVS)  Large  Amplitude  Simulator/Wide  Angle  Visual  System.  The 
comparison  is  conducted  for  the  50,000  pounds  gross  weight  to 
establish  the  degree  of  validity  associated  with  the  pre¬ 
selected  analytical  approaches. 

The  analytical  results  indicate  that  the  tandem  configurations 
(up  to  200,000  pounds)  can  hold  position  with  essentially  the 
same  accuracy  regardless  of  the  specific  gross  weight.  This 
property  is  the  result  of  the  overwhelming  dominance  of  the 
automatic  control  over  the  inherent  aircraft  characteristics, 
as  well  as  of  the  basic  dynamic  similarity  of  all  gross  weight 
configurations . 

Analysis  of  the  hover  hold  capability  in  turbulence  reveals 
that  the  piloted  performance  is  mainly  a  function  of  specific 
control  augmentation  modes  (command  systems  or  control  types) , 
and  is  only  weakly  dependent  on  actual  gross  weight  configura¬ 
tion.  This  is  another  indication  of  the  dominance  of  the 
automatic  pilot  assist  modes  over  the  inherent  dynamic  response 
characteristics  of  helicopters  in  hover.  In  general,  a  good, 
precise,  hover  hold  capability  is  obtainable  with  horizontal 
velocity  command  system  involving  very  little  attitude  re¬ 
sponse. 

The  simulator  evaluation  of  the  50,000-pound  configuration 
confirms,  in  general,  the  basic  analytical  results  relating  to 
the  piloted  and  the  automatic  system  performance.  In  case  of 
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piloted  performance,  the  analysis  predicts  well  the  best  ob¬ 
tainable  hover  hold  capability  for  all  command  systems,  but 
the  average  actual  performance  on  simulator  agrees  exactly 
with  the  analytical  predictions  only  for  linear  velocity  type 
command  systems.  The  difficulty  in  predicting  piloted  per¬ 
formance  with  marginal  command  systems  (e.g. ,  attitude  or  ac¬ 
celeration  types)  stems  from  the  uncertainty  in  assessing  the 
pilot  apprehension  level  resulting  from  display  limitations, 
damping  level  preferences,  and  indirect  visual  and/or  physical 
cues  (e.g., attitude  changes  during  position  control;  rapid 
position  display  motion  because  of  small  scale,  etc). 

The  general  results  of  this  research  effort  provide  a  vali¬ 
dated  basis  for  further  analysis  and  piloted  simulation  of 
VHLH  configurations  regarding  their  performance  capabilities, 
handling  qualities  criteria  and  automatic  system  requirements 
in  hover  and  forward  flight.  The  package  containing  MLH  and 
VHLH  configuration  data  provides  the  aircraft  data  for  analyt¬ 
ical  studies.  This  data  package,  combined  with  the  validated 
simulator  model  at  the  Northrop  Corporation,  Aircraft  Division, 
LAS/WAVS  facility,  is  suitable  for  further  piloted  simulator 
evaluation  of  hover  and  forward  flight  helicopter  performance 
with  and  without  externally  slung  loads.  Finally,  the  analyt¬ 
ical  approach  used  here  and  validated  against  the  piloted 
simulator  results  enables  the  designer  to  optimize  his  low 
speed  control  concepts  inexpensively  with  a  significant  con¬ 
fidence  level. 
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FOREWORD 


This  report  constitutes  the  work  of  the  Boeing  Vertol  Company 
in  performance  of  U.S.  Army  Air  Mobility  Research  and  Develop¬ 
ment  Laboratory  Contract  DAAJ02-72-C-0095  (DA  Task  1F162204- 
AA4401) .  The  work  was  performed  from  June  1972  through  August 
1973  with  Mr.  James  M.  Davis  as  the  Program  Manager,  and 
Messrs.  R.T.  Lytwyn  and  F.  White  as  the  principal  investiga¬ 
tors.  The  Army  technical  representative  was  Mr.  Robert  P. 
Smith  of  Eustis  Directorate,  USAAMRDL,  Fort  Eustis,  Virginia. 

The  ground-based  simulator  evaluation  work  was  performed  at 
the  Northrop  Corporation,  Aircraft  Division  Aerosciences 
Laboratory  in  Hawthorne,  California.  The  simulator  mechani¬ 
zation  was  accomplished  under  the  guidance  of  John  B.  Sinacori. 
Mr.  A.  Murakoshi  was  the  principal  application  engineer,  with 
Messrs.  William  Ross  and  Raymond  Silvestri  assisting  him 
throughout  the  program. 

Mr.  A. P .  Santa-Maria  from  the  Boeing  Vertol  Experimental 
Flight  Operations  was  the  simulator  test  pilot  on  this  pro¬ 
gram,  and  Mr.  E.  Low  from  the  Aerodynamics  Technology  Staff 
performed  the  configuration  definition  work. 
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INTRODUCTION 


BACKGROUND 


Potential  applications  for  large  helicopters  with  external 
load  capability  include  the  stacking  of  containers  on  trucks 
or  ships,  bridge  placement,  and  modular  housing  construction. 

A  high  degree  of  precision  is  required  from  the  helicopter  if 
these  types  of  operations  are  to  be  performed  in  a  timely  and 
safe  manner. 

The  ability  of  a  helicopter  to  hover  accurately  with  respect 
to  a  fixed  reference  point  depends  on  the  inertial  and  aero¬ 
dynamic  characteristics  of  the  aircraft,  the  design  of  the 
automatic  flight  control  system,  and  the  response  of  the  rotor 
and  flight  control  system.  At  present,  no  experience  with 
large  helicopters  J( 80 ,000  pounds  and  larger  gross  weights) 
exists  which  would  indicate  their  potential  capabilities  or 
limitations  during  precise  hover  hold  missions.  As  the  size 
of  the  helicopter  is  increased,  larger  rotors  are  required 
which  must  operate  at  lower  rpm.  This  leads  to  lower,  rotor- 
induced  frequencies,  which  might  conceivably  alter  the  basic 
responses  normally  experienced  with  existing  medium-gross- 
weight  aircraft.  Pure  extrapolation  to  a  larger  size  contains 
a  sufficient  number  of  unknown  factors  affecting  the  flying 
qualities  to  warrant  a  closer  analytical  and  simulator  evalua¬ 
tion  of  the  performance  of  such  aircraft  during  precise  hover 
hold  missions.  The  questions  of  aircraft  size  effect  on  the 
pilot's  ability  to  perform  these  missions  must  be  addressed 
in  analyses  and  simulations.  The  control  system  pilot-assist 
requirements,  as  well  as  the  automatic  hover  hold  functions, 
must  be  synthesized  and  evaluated. 

The  sensitivity  of  these  helicopters  to  turbulence  levels  must 
also  be  established  if  an  all-weather  operational  capability 
is  to  be  ensured.  The  acceleration  environment  at  the  crew 
station  should  be  investigated  since  the  cockpit  position  with 
respect  to  aircraft  center  of  gravity  will  increase  with  the 
helicopter  size,  and  too  large  acceleration  levels  will  inter¬ 
fere  with  the  precise  piloting  task. 

OBJECTIVES  AND  SCOPE 

The  main  objective  of  the  analytical  portion  of  this  investi¬ 
gation  is  the  evaluation  of  helicopter  size  effects  on  its 
capability  to  perform  precise  hover  position  hold  missions. 

For  this  purpose,  the  study  defines  parametrically  four 
tandem-rotor  helicopter  configurations  with  the  gross  weight 
capability  of  50,000,  80,000,  120,000,  and  200,000  pounds, 
and  synthesizes  the  automatic  control  system  functions  for 
precise  hover  capability  and  pilot  control  augmentation  (as¬ 
sist  functions) .  A  single  precision  hover  system  (PHS)  and 
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stability  augmentation  system  (SAS)  is  used  for  all  gross 
weight  configurations.  Several  pilot-assist  functions  are 
designed  using  different  combinations  of  longitudinal  cyclic 
pitch  (DCP) ,  with  conventional  lateral  cyclic  control.  Auto¬ 
matic  feedback  functions  are  also  included  to  provide  the 
pilot  with  direct  acceleration,  velocity,  or  altitude  control 
on  the  stick.  A  full-flight  envelope  mathematical  model  is 
utilized  to  perform  piloted  evaluation  and  to  compare  the 
results  with  analytical  predictions.  The  hover  position  hold 
capabilities  and  pilot  station  accelerations  are  defined  as  a 
function  of  turbulence  level,  pilot  control  augmentation  modes, 
and  gross  weight  configurations.  The  unpiloted  aircraft  posi¬ 
tion  hold  capabilities  with  PHS  augmentation  are  compared 
with  the  piloted  performance  without  ground  position  and  ve¬ 
locity  hold.  Several  test  points  are  compared  directly  with 
the  piloted  simulator  results  in  terms  of  turbulence  effects, 
control  assist  functions,  and  pilot  station  acceleration 
levels . 
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DEVELOPMENT  OF  PARAMETRIC  CHARACTERISTICS  FOR 
ADVANCED  TANDEM- ROTOR  HEAVY  LIFT 


HELICOPTER  CONFIGURATIONS 


The  four  tandem-rotor  helicopters  defined  in  this  section  are 
based  on  scaling  the  existing  HLH  designs  up  into  the  VHLH 
gross  weight  region.  The  ground  rules  for  scaling  of  the 
aerodynamic  performance  and  fuselage  aerodynamic  forces  and 
moments  are  discussed  below,  and  the  dynamic  characteristics 
of  the  basic  aircraft  are  presented  on  page  13. 

ADVANCED  TANDEM-ROTOR  HEAVY  LIFT  HELICOPTER  (HLH) 
CONFIGURATIONS 


Sizing  Ground  Rules  and  Aircraft  Parameters 

The  four  HLH  configurations  shown  in  Figure  1  are  based  on  the 
scaling  of  the  Boeing-Vertol  Model  301  Heavy  Lift  Helicopter, 
with  a  design  gross  weight  of  from  120,000  pounds  to  the 
equivalent  of  tandem-rotor  helicopters  at  higher  and  lower 
gross  weights.  As  such,  they  are  based  on  the  following  scal¬ 
ing  ground  rules: 

1.  Disc  loading  =  8.9  lb/ft 

2.  Rotor  overlap  =  23.8%  of  rotor  diameter 

3.  Rotor  solidity  =  0.0923 

4.  Tip  speed  =  750  fps 

5.  All  reference  lengths,  except  the  main  landing  gear, 
are  proportional  to  design  gross  weight.  The  main 
landing  gear  is  assumed  to  be  constant  in  order  to 
accommodate  an  8-by-8-by  20-ft  cargo  container. 

6.  All  reference  areas,  except  the  equivalent  flat  plate 
area  (fe),  are  proportional  to  design  gross  weight. 

The  flat  plate  area,  (fe),is  based  on  the  following: 

•  Main  landing  gear  flat  plate  area  varies  as  the 
square  root  of  gross  weight. 

•  All  other  items  vary  linearly  with  gross  weight. 

7.  Engine  and  transmission  ratings  are  proportional  to 
design  gross  weight  . 

8.  Moments  of  inertia  are  scaled  from  the  Model  301  Heavy 
Lift  Helicopter  data,  using  the  length  and  area  pro¬ 
portionality  described  above  . 
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NORMAL  GROSS  WT: 


50,000  LB 
R  =  29.90  FT 
Dr  =  45.57  FT 
h  =  6.0  FT 
hH  =  6.83  FT 


80,000  LB 
R  =  37.82  FT 
Dr  =  57.64  FT 
h  =  7.6  FT 
hH  =  8.63  FT 


Figure  1.  Advanced  Tandem-Rotor 
HLH  Configuration  Data. 
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The  aircraft  parameters  defining  each  design  gross  weight 
configuration  are  summarized  in  Table  I. 

The  disc  loading  of  8.9  pounds  per  square  foot  represents  more 
realistically  the  current  tandem-rotor  helicopter  designs, 
while  attempting  to  keep  the  downwash  velocity  sufficiently 
small  at  reasonable  power  level  requirements. 

Fuselage  Aerodynamics  and  Scaling  Factors 


Tables  II  through  VII  summarize  the  fuselage  aerodynamic 
forces  and  moments  for  the  baseline  HLH  configuration  at 
120,000  pounds  gross  weight.  The  tables  are  based  on  wind 
tunnel  test  results  for  a  solid  model  of  a  tandem-rotor  heavy 
lift  helicopter.  Data  presented  in  the  tables  are  functions 
of  fuselage  angle  of  attack,  a,  and  fuselage  sideslip  angle, 

6.  Tabular  force  and  moment  information  is  expressed  in  the 
body  axis  in  the  form:  drag  (D/qd) ,  sideforce  (Y/qd) ,  lift 
(L/qd)  ,  pitching  moment  (M/qd) ,  rolling  moment  (^f/qd)  ,  and  yawing 
moment  (N/qd) ,  with  a  from  -90°  to  +90°,  and  8  from  -90°  to 
+90°.  From  90°  to  180°,  the  curves  are  assumed  to  be  mirror 
image  values  of  0°  to  90°. 

Moment  data  is  expressed  about  a  center-of-gravity  location  at 
which  the  wind  tunnel  model  was  tested.  Moment  corrections 
proportional  to  the  distance  to  the  c.g.  location  being  flown 
with  the  full  force  simulation  model  are  applied  to  the  stored 
data. 


Tabular  fuselage  data  usually  does  not  account  for  rotor  hub 
drag,  momentum  drag,  and  leakage  drag  and  therefore  must  be 
adjusted  to  the  correct  level  before  being  used.  A  Af  cor¬ 
rection  is  applied  to  the  tabulated  lift,  drag,  and  sice  force 
data  to  account  for  these  items. 


The  fuselage  aerodynamic  data  shown  in  Tables  II  through  VII 
are  already  corrected  for  the  rotor  hub  drag,  and  are  based 
on  the  total  equivalent  flat  plate  area  of  138.0  sq  ft.  This 
represents  very  closely  the  120,000-pound  configuration,  and 
the  tables  can  be  used  directly  for  that  configuration  gross 
weight.  For  configurations  other  than  the  120,000  pounds,  the 
following  scale  factors  should  be  applied  to  arrive  at  the 
corresponding  fuselage  aerodynamic  tables: 


qdl 


TABLE  II 


TABLE  IV 
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TABLE  I.  ADVANCED  TANDEM-ROTOR  HLH  CONFIGURATION  PARAMETERS 
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TABLE  II.  BASELINE  TANDEM- ROTOR  HLH  FUSELAGE  DRAG  TABLES,  D/qd 
REFERENCE:  STA.  596",  WATER  LINE  164" 

120,000-POUND  GROSS  WEIGHT 
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TABLE  IV.  BASELINE  TANDEM-ROTOR  HLH  FUSELAGE 
SIDEFORCE  TABLES,  Y/qd 
120,000-POUND  GROSS  WEIGHT 
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TABLE  V.  BASELINE  TANDEM-ROTOR  HLH  FUSELAGE  PITCH 
MOMENT  TABLES,  M/qd 
120,000-POUND  GROSS  WEIGHT 
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TABLE  VI.  BASELINE  TANDEM- ROTOR  HLH  FUSELAGE  ROLL 
MOMENT  TABLES, -t/qd 
120,000-POUND  GROSS  WEIGHT 
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TABLE  VII.  BASELINE  TANDEM-ROTOR  HLH  FUSELAGE 
YAW  MOMENT  TABLES,  N/qd 
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L/9d  (fe)  (L/q<J)TABLE  III 
M/qd  ‘ik^)  (M/qd)  table  V 
%d  =(  ^ref)  (  /qd) table  VI 
N/qd  =  (^7)  (N/qd)  TABLE  VII 

where 

£eREP  =138.0  ft2 

f  =  tabulated  in  TABLE  I 

e 

Rref  -  45.0  ft 

R  =  tabulated  in  TABLE  I 


These  scaling  formulas  are  based  on  the  assumptions  that  all 
non-dimensional  aerodynamic  force  and  moment  coefficients  are 
constant  for  each  gross  weight  configuration.  Only  the  scal¬ 
ings  of  reference  lengths  and  areas  are  required,  and  these 
are  proportional  to  rotor  radius  and  (rotor  radius) ^  respec¬ 
tively. 

DYNAMIC  CHARACTERISTICS 

The  stability  derivations  presented  in  this  section  are  based 
on  the  standard  Boeing-Vertol  trim  and  stability  digital  pro¬ 
gram  A-97 .  The  program  computes  the  helicopter  trim  by 
means  of  an  iterative  solution  to  the  six  steady-state  equa¬ 
tions,  summing  all  forces  and  moments  about  the  fixed  body 
axis.  A  numerical  approach  to  the  solution  of  the  rotor  flap¬ 
ping  and  force  equations  is  used,  allowing  for  the  automatic 
accounting  of  blade  stall,  reverse  flow,  and  compressibility 
effects.  The  simplifications  or  assumptions  contained  in  the 
rotor  analysis  include  the  following: 

•  Induced  velocity  distribution  is  assumed  to  be  uniform. 

•  Blade  lag  and  all  elastic  degrees  of  freedom  are 
neglected. 
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•  Nonsteady  aerodynamic  and  spanwise  flow  effects  are 
neglected. 

The  stability  derivatives  and  control  powers  are  then  obtained 
numerically  by  evaluating  the  changes  in  the  helicopter  forces 
and  moments  from  their  trim  value  caused  by  a  small  perturba¬ 
tion  in  an  independent  parameter.  These  derivatives  are 
quasi-static;  i.e.,  the  rotor  is  assumed  to  reach  its  new 
steady-state  condition  instantaneously  as  the  perturbation  is 
applied. 

The  stability  derivatives  data  are  presented  for  sea-level 
standard  conditions  (midpoint  c.g.  position)  in  hover  at  all 
design  gross  weights  and  all  estimated  empty  gross  weights. 

An  analysis  is  then  made  of  the  effects  of  scaling  on  impor¬ 
tant  hover  derivatives  and  the  associated  dynamic  properties. 
All  control  sensitivity  data  are  based  on  constant  cockpit 
control-to-blade  pitch  angle  mechanical  ratios  consistent  with 
the  typical  tandem-rotor  HLH  configuration. 

Stability  Derivatives 

A  summary  of  stability  derivatives  for  all  design  gross  weight 
configurations  in  hover  is  shown  in  Table  VIII.  The  most  im¬ 
portant  derivatives,  which  define  the  helicopter  dynamics  in 
hover,  are  plotted  in  Figures  2  and  3.  The  following  discus¬ 
sion  elucidates  the  scaling  properties  of  the  major  derivatives 
plotted  in  these  figures. 

Longitudinal  Axis 

The  most  important  derivatives  in  hover  are  represented  by  the 
vertical  damping,  Zw,  the  pitch  damping,  Mq,  and  the  pitch 
acceleration  due  to  longitudinal  velocity,  Mu.  These  deriva¬ 
tives,  together  with  the  control  sensitivities  Z5 c  and  M<sB, 
determine  the  aircraft  response  to  pilot  input  in  hover. 

The  vertical  damping ^  Zw»  is  determined  by  the  sensitivity  of 
rotor  thrust  to  vertical  velocity.  According  to  Reference  1, 
this  derivative  in  hover  can  be  approximated  by  the  formula 
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TABLE  VIII.  ADVANCED  TANDEM- ROTOR  HLH 
STABILITY  DERIVATIVES 

Hover,  midpoint  C.G.,  sea  level,  standard  day 
7S0  ft/sec  tip  speed 


Param. 

Symbol 

Units 

9o 

deg 

bitf 

deg 

bitr 

deg 

Xu 

"sec 

xw 

"sec 

Xq 

f"sec 

X«B 

ft/sec2 

X«c 

msm  -is 

BHEv 

in. 

Zu 

"sec 

•Zw 

"sec 

Zq 

f"sec 

Z«B 

f"sec2 

50,000  Lb 


Configuration 

80,000  Lb  120,000  Lb 


^sec1 

in. 

^sec1 


"sec2- 


.2926 

.01010 


Bnpty 

Gross 

Weight 

Design 

Gross 

Weight 

Empty 

Gross 

Weight 

2.804 

2.805 

2.769 

-.40 

-.40 

-.40 

-2.80 

-2.80 

1 

to 

CO 

o 

-.01260 

-.01756 

-.01180 

.01389 

.00794 

.01476 

-.06072 

.7521 

-.2741 

.2303 

.1-91 

.  2448 

.4631 

.  2925 

.4950 

.01762 

.01013 

.01876 

-.3556 

-.2229 

-.3771 

-1.5102 

-1.2620 

-2.0971 

.1418 

.1196 

.1394 

-11.478 

-7 . 308 

-12.151 

.00457 

.00396 

.00361 

-.00401 

-.00335 

-.00314 

-.8166 

-.8053 

-.8400 

.  2946 

.2241 

.2401 

-.1213 

-.1031 

-.0949 

-.01472 

-.01892 

-.01395 

-.  4400 

-1.1161 

-.4147 

-.01400 

-.03977 

-.00541 

1.  3372 

1.3143 

1.3164 

-.08461 

-.05521 

-.07885 

-.00618 

-.00672 

-.00526 

-.6473 

-.6508 

-.7676 

.02171 

.036  31 

.01354 

.6412 

.4813 

.  5838 

-.209] 

-.1852 

-.1891 

.0000271 

.0000556 

0000115  . 

-.001963 

b. 000444 

-.00371 

-.04756 

-.08608 

-.04190 

.000905 

.000785 

.000966 

.1541 

.1796 

.1116 

1920  .2 


200,000  Lb 


Design  Empty 
Gross  Gross 
Weight  Weight 


-0.40 

-2.80 

-.01088 

.01577 

-.■>■>51 


.1440  .2655 


.3015 

.00995 

-.2263 

-1.9504 


.  S341 
.01815 
-.3975 
-3.295 


-7.397  -12.881 

.00266  .00243 

-.00217  -.00196 

-.8900  -.9209 

.1561  .1642 


-.0646 

-.01827 

-1.6800 

-.0380 


-.0572 

-.01335 

-.3962 

.0238 


1.2964  1.2942 

-.05830  -.08427 

-.00535  -.00400 


-.9137 

.05379 


-1.1346 

.01423 


.3931  .4788 

-.1608  -.1531 

.0000272  .0000154 
.000309  -.00397 

3903  -.08443  -.03858 

.000435  .000719 

.1119  .0657 


Denotes  the  most  important  stability  derivatives 


50,000  80,000  120,000  200,000 


DESIGN  GROSS 
WEIGHT 

WEIGHT  EMPTY 
(LOW  DISC 
LOADING) 


CONFIGURATION  GROSS  WEIGHT  -  LB 


Figure  2.  Major  Longitudinal  Stability 
Derivatives  in  Hover. 
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DESIGN  GROSS 
WEIGHT 

WEIGHT  EMPTY 
(LOW  DISC 
LOADING) 


Figure  3.  Major  Lateral-Directional  Stability 
Derivatives  in  Hover. 
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Since  the  tip  speed,  \'T,  the  disc  loading,  and  the  blade-chord 
to-rotor  radius  ratio  are  constant  with  the  configuration 
gross  weight,  the  vertical  damping  remains  nearly  constant 
with  gross  weight  (see  Figure  2) .  The  factor 


Disc  Loading 
|  NPa  VT2 


so  at  lower  disc  loading,  the  vertical  damping  will  increase 
for  a  given  configuration.  This  is  really  the  reflection  of 
decreasing  helicopter  mass  as  the  gross  weight  of  a  given 
configuration  is  reduced.  The  slight  increase  in  Zw  for  the 
empty  weight  configuration  with  higher  design  gross  weights 
represents  the  smaller  ratio  of  empty-to-design  gross  weight 
factor  for  heavier  helicopters. 

The  pitch  damping,  Mg,  in  hover  is  proportional  to  the  verti¬ 
cal  damping,  so 

Mq  =  ZW  /ija  ) 

\  yy/ 

Again,  since  Zw  is  constant  with  configuration  design  gross 
weight,  Mq  should  scale  as  the  factor  in  brackets.  For  the 
particular  inertia  and  aircraft  dimension  scaling  noted  in 
Table  I,  the  scaling  factor  increases  by  about  15%  from  the 
50,000-pound  to  the  200,000-pound  configuration.  This  scaling 
property  is  reflected  in  the  increased  Mq  with  configuration 
design  gross  weight.  At  the  empty  gross  weight  (low  disc 
loading) ,  the  increased  Mq  is  due  to  lower  pitch  inertia. 

The  ipitch  acceleration  due  to  forward  speed,  Mu,  can  be  ap¬ 
proximated  in  hover  (Reference  1)  by  the  formula 


so  Mu  will  scale  according  to  the  factor  (herm/Iyy) .  The  air¬ 
craft  dimensional  and  inertia  scaling  shown  in  Table  I  results 
in  a  decreasing  Mu  with  the  configuration  gross  weight.  The 
slight  decrease  in  Mu  at  empty  gross  weight  is  due  to  lower 
disc  loading  and  lower  pitch  inertia. 
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The  vertical  control  sensitivity,  Z^r?/  is  proportional  to  the 
rotor  thrust  derivative  with  respect  to  the  blade  pitch  set- 


_0)  ('-  r  rj 

so  the  vertical  control  sensitivity,  Z$c,  remains  nearly  con¬ 
stant  for  all  design  gross  weight  configurations  (Figure  2) 
and  increases  with  the  lower  disc  loading  (empty  gross  weight) 
The  increased  c  at  low  gross  weight  is  again  the  direct  re¬ 
sult  of  lower  mass  at  constant  thrust  derivative  with  respect 
to  collective  pitch. 

The  longitudinal  stick  sensitivity,  M^b,  is  directly  propor- 
tional  to  Z6C  in  tandem-rotor  helicopters,  because  the  pitch 
control  in  tandem-rotor  helicopters  is  derived  mainly  from 
differential  change  in  blade  collective  pitch  angle  on  the 
forward  and  rear  rotors.  The  variation  of  M6b  with  aircraft 
size  can  then  be  expressed  in  the  proportional  form 

m«b  -  zec 

and  Msg  will  scale  as  the  factor  (£m/Iyy) .  For  the  particular 
parameters  shown  in  Table  I,  will  decrease  by  more  than 
40%  as  the  aircraft  configuration  gross  weight  is  increased 
from  50,000  pounds  to  200,000  pounds. 

Lateral-Directional  Axis 


ting.  In  hover  (Reference  1)  , 


z«o  “  z«c  s  -  3 


The  most  important  lateral-directional  derivatives  in  hover 
are  represented  by  the  roll  damping,  Lp,  the  yaw  damping,  Lr, 
and  the  roll  acceleration  due  to  side  velocity,  Ly  These 
derivatives,  together  with  the  control  sensitivities,  Lgs  and 
N^r,  define  the  basic  aircraft  response  to  the  pilot  input  in 
hover. 

The  roll  damping,  Lp ,  is  a  function  of  rolling  moment  due  to 
rotor  disc  lateral  tilting  brought  about  by  aircraft  roll  rate 
An  approximate  expression,  valid  in  hover,  is  given  by 


^erWG  + 


xx 


V„ 


2 


As  shown  in  Table  I,  VT  is  constant  and  Y  is  nearly  constant 
for  all  gross  weight  configurations.  Lp  will  then  scale  as 
the  rotor  radius  multiplied  by  the  expression  in  brackets. 

The  expression  in  brackets  is  plotted  in  Figure  4  in  terms  of 
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THRUST  TILT  CONTRIBUTION 


T ,  her  Wq 
laic  =  ixx 


CENTRIFUGAL  MOMENT  CONTRIBUTION ; 


T  ..  _  NoBVt2  (  e3\ 

laiC-rixx  \r) 


Wfe  =  TOTAL  GROSS  WEIGHT 


-  DESIGN  GROSS  WT 

-  EMPTY  GROSS  WT 
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Figure  4.  Average  Roll  Acceleration  Sensitivities 
to  Rotor  Disc  Tilt  in  Hover. 
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the  rotor  thrust  tilt  contribution,  an<*  the  centrifugal 

moment  part,  L^ic,  an<*  as  the  sum  (Laic  +  LAlc)  •  the  total 
roll  acceleration  sensitivity.  As  shown  in  the  figure,  the 
total  sensitivity  decreases  with  the  configuration  gross 
weight,  but  the  decrease  is  much  less  than  the  inverse  ratio 
of  rotor  radii.  Lp  will  then  increase  with  configuration 
gross  weight  as  shown  in  Figure  3.  At  empty  gross  weights, 
the  higher  Lp  is  a  reflection  of  the  dominant  contribution  of 
the  centrifugal  forces  to  the  rolling  moment,  which  scales  as 
the  inverse  ratio  of  roll  inertia. 

The  yaw  damping,  Nr,  is  related  to  the  rotor  disc  tilt  brought 
about  by  lateral  translational  velocities  at  forward  and  aft 
rotors.  As  such,  it  is  a  direct  function  of  Yv  of  each  rotor, 
which  can  be  approximated  in  hover  by  the  formula 


At  constant  disc  loading  and  constant  (C/R) ,  Yv  will  remain 
constant  for  all  gross  weight  configurations.  At  lower  disc 
loading  (e.g.,  weight  empty),  Yv  will  decrease  proportionately 
to  the  disc  loading.  The  yaw  damping  will  then  scale  with 
configuration  gross  weight  as  the  ratio  U2 m/Izz)  ,  which,  for 
the  parametric  aircraft  definition  given  in  Table  I,  is  equal 
to  unity  at  the  design  gross  weight,  and  somewhat  less  than 
unity  at  the  empty  gross  weights.  Consequently,  the  yaw  damp¬ 
ing,  Nr,  remains  nearly  constant  with  the  design  gross  weight 
and  is  decreased  uniformly  at  lower  disc  loading. 

The  roll  acceleration  due  to  sideward  velocity,  Ly,  is  also 
directly  related  to  the  rotor  disc  tilt  brought  about  by  the 
side  velocity.  As  such,  its  dependence  on  design  gross  weight 
configurations  can  be  related  to  Yv  in  the  proportional  form 


The  scaling  factor  (h^rm/Ixx)  decreases  with  the  design  gross 
weight,  while  Yv  remains  constant.  The  derivative  Lv  will 
then  decrease  at  higher  design  gross  weights.  It  will  also 
decrease  proportionately  with  the  disc  loading. 


21 


The  lateral  cyclic  control  sensitivity,  L^c«  is  a  direct  func- 
tion  of  the  lateral  rotor  disc  tilt  per  unit  of  cyclic  pitch 
angle.  In  general,  the  rolling  moment  per  unit  cyclic  pitch 
comprises  two  most  powerful  contributions:  the  lateral  com¬ 
ponent  of  the  thrust,  and  the  hub  moment  due  to  the  offset 
flapping  hinges.  At  hover,  L6S  can  be  approximately  defined 
in  terms  of  roll  moment  per  unit  lateral  cyclic  in  the  form 


L*s 


AIC 


The  first  term  in  the  brackets  represents  the  contribution  of 
the  thrust  tilt,  while  the  second  is  the  centrifugal  force 
moment.  The  first  factor  will  scale  with  the  design  gross 
weight  as  (h?rWG/Ixx) ,  and  will  therefore  decrease  with  in¬ 
creasing  design  gross  weight,  diminishing  slightly  at  empty 
gross  weights  (see  Figure  4) .  The  second  factor  remains  con¬ 
stant  at  all  configuration  gross  weights,  but  scales  as  (  ^ ^ 

V*xx  J 

from  normal  gross  weight  to  weight  empty.  This  centrifugal 
force  contribution  will  remain  essentially  constant  with  disc 
loading. 

The  yaw  control  sensitivity,  N6r,  for  a  tandem-rotor  helicopter 
is  generated  by  the  lateral  cyclic  pitch  applied  differentially 
on  forward  and  aft  rotors.  As  such,  it  can  be  expressed  by 
the  following  approximate  formula,  based  on  thrust  tilt,  above: 


N6r  «  Nmc  =  *WG/IZZ 


N<$r  will  then  scale  as  the  factor  (^Wq/Izz)  ,  and  will  there¬ 
fore  decrease  with  the  increased  configuration  gross  weight. 
The  diminishing  property  of  N6r  at  empty  gross  weight  (see 
Figure  3)  is  attributable  to  the  particular  geometric  and 
inertial  scaling  properties  of  the  configurations,  rather  than 
to  disc  loading  effects. 

Dynamic  Properties  of  Basic  Aircraft  Configurations 

The  dynamic  characteristics  discussed  in  this  section  are 
based  on  aircraft  rigid-body  behavior  as  defined  by  the  un¬ 
coupled  longitudinal  and  lateral-directional  stability  de¬ 
rivatives  listed  in  Table  VIII.  An  approximation  of  the 
characteristic  roots  in  hover,  based  on  Reference  1,  is  used 
to  elucidate  the  variations  in  these  roots  with  configuration 
gross  weight  and  disc  loading.  The  complete  dynamic  charac¬ 
teristics,  without  simplifications,  are  tabulated  in  Table  IX. 
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TABLE  IX.  LATERAL- DIRECTIONAL  AND  LONGITUDINAL 
ROOTS  IN  HOVER 
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Longitudinal  Dynamics 

The  longitudinal  dynamics  in  hover  can  be  approximated  by  the 
following  expressions  for  the  corresponding  characteristic 
roots: 

Vertical  Mode:  rw  =  Zw 

Pitch  Subsidence  Mode:  rq  =  Mq 

Pitch  -  Longitudinal  Oscillations  (Mq  <  0) : 


i 


Vertical  Mode :  This  mode  determines  the  aircraft  vertical 
velocity  response  to  collective  stick,  with  the  time  constant 
equal  to 


Figure  2  shows  that  the  vertical  response  at  normal  gross 
weight  remains  constant  for  all  configurations  with  a  time 
constant  of  4.5  seconds.  At  weight  empty,  the  vertical  time 
constant  will  shorten  to  less  than  3.0  seconds.  Such  time  con¬ 
stants  are  too  long  for  precise  vertical  control,  and  collec¬ 
tive  control  augmentation  will  be  required  for  precise  hover 
tasks. 

Pitch  Subsidence  Mode:  As  shown  in  Table  IX,  the  pitch  subsi- 
dence  mode  remains  nearly  constant  for  all  gross  weight  con¬ 
figurations  with  a  time  constant  of  1  second.  Pitch  damping 
augmentation  will  be  required  in  order  to  improve  pure  pitch 
control  and  reduce  the  pitch  time  constant  to  a  reasonable 
magnitude. 

Longitudinal  Oscillatory  Mode:  This  mode  does  not  determine 
the  actual  aircraft  response  to  pilot  inputs,  but  its  unstable 
character  is  generally  detrimental  to  control  in  any  turbu¬ 
lence  level.  The  frequency  and  the  damping  ratio  of  this 
mode  are  functions  of  the  ratio  of  Mu  to  Mg  alone,  so  the  mode 
changes  very  little  with  gross  weight  configuration  and  disc 
loading.  The  unstable  feature  of  this  mode  in  all  hovering 
helicopters  is  easily  suppressed  by  pitch  damping  augmentation 
through  the  differential  collective  pitch  feedback. 
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Lateral-Directional  Dynamics 


The  lateral-directional  dynamics  in  hover  can  be  described  by 
the  following  approximate  expressions  for  the  characteristic 
roots  (poles) : 


Roll  Subsidence  Mode:  r _ 

rs 


Spiral  (Yaw)  Mode:  r  =  N 

sp  r 


Dutch  Roll  Mode: 


+ 


i 


The  roll  subsidence  mode  is  the  primary  response  mode  around 
which  the  pilot  controls  the  helicopter  in  roll  and  sideward 
motion.  As  shown  in  Table  IX,  the  time  constant  of  this  mode, 
equal  to  | l/rrs | ,  remains  between  1.1  and  0.9  second  for  all 
gross  weight  configurations.  As  such,  it  is  generally  too 
long  for  precise  control,  so  roll  rate  feedback  is  generally 
required  to  improve  the  roll  control  on  tandem,  fully  articu¬ 
lated  rotors. 


The  spiral  or  yaw  mode  in  hover  determines  the  aircraft  re¬ 
sponse  to  pedals.  In  tandem-rotor  helicopters  this  mode  is 
nearly  neutrally  stable  (see  Table  IX)  because  of  insignificant 
directional  stability  in  hover.  For  precise  control,  all  gross 
weight  configurations  will  require  stability  augmentation  in 
the  yaw  axis . 

The  Dutch  roll  mode  is  an  unstable  oscillatory  mode  dominated 
by  roll  and  sideward  motion  with  some  yaw.  The  unstable  char¬ 
acter  of  this  motion  is  generally  detrimental  to  aircraft  con¬ 
trol  in  any  turbulence  in  hover.  In  hover  the  mode  can  be 
stabilized  usually  by  a  tight  roll  feedback  loop  (rate  and 
attitude) . 
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FULL-FLIGHT  ENVELOPE  MATHEMATICAL  MODEL 


The  airframe  and  the  control  system  analytical  models  outlined 
in  this  section  are  the  same  as  the  standard  Boeing-Vertol 
model  utilized  for  the  piloted  simulation  of  the  Model  347 
helicopter  and  for  Heavy  Lift  Helicopters.  As  such,  the  model 
contains  all  the  basic  features  of  tandem-rotor  helicopters, 
including  a  simplified  external  load  suspension  system  utiliz¬ 
ing  two  longitudinally  disposed  cables.  Only  the  basic  fea¬ 
tures  of  the  model  are  described  in  this  report.  A  detailed 
description  of  the  model  is  contained  in  Reference  5. 

AIRFRAME  DYNAMIC  MODEL 


The  airframe  and  rotor  dynamic  model  used  for  the  full-flight 
envelope  tandem  helicopter  simulation  is  described  in  this 
section.  To  illustrate  how  this  model  works,  brief  descrip¬ 
tions  of  each  major  calculation  have  been  prepared  and  are 
presented  below  in  the  general  order  in  which  they  are  per¬ 
formed,  starting  with  equations  of  motion  and  working  back 
around  to  the  forces  and  moments  which  drive  these  motion 
equations . 

Figure  5  illustrates,  in  block  diagram  format,  how  the  various 
calculations  in  the  model  are  related.  Shown  at  the  far 
right-hand  side  of  the  chart  are  the  equations  of  motion  being 
fed  by  the  summed  forces  and  moments.  Other  portions  of  the 
diagram  are  associated  with  the  (  omputational  steps  described 
be  low . 

Airframe  Calculation  Summary 

1.  Basic  equations  of  motion  are  solved,  providing  linear 
and  angular  accelerations  about  all  three  axes.  Euler 
angle  rotations  are  performed  to  orient  the  airframe 
in  inertial  space. 

2.  Remote  velocity  is  resolved  into  the  body  axis  to  pro¬ 
vide  fuselage  angle  of  attack  and  sideslip  information. 
Effect  of  wind  on  velocity  resolution  is  considered, 

as  is  forward  rotor  downwash  interference  on  fuselage. 

3.  Remote  velocity  is  resolved  through  rotor  shaft  inci¬ 
dence  angle  (and  is  corrected  for  aircraft  pitching 
motion)  to  derive  rotor  shaft  normal  plane  (S.N.P.) 
wind  axis  velocities  and  sideslip  (g)  values. 

4.  Rotor  advance  (P )  and  inflow  (X)  ratios  are  computed. 
Inflow  ratio  calculations  account  for  rotor-on-rotor 
interference. 
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Figure  S.  Tandem-Rotor  Helicopter  Aircraft 
Mathematical  Model. 
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aft 


5.  Control  axis  angle  of  attack  (acA)  for  each  rotor  is 
calculated.  Rotor-rotor  interference  effects  are 
included. 

6.  Rotor  input  collective  and  cyclic  controls  (0O,  B iq, 
Aic)  are  compiled,  and  both  longitudinal  and  lateral 
cyclics  are  resolved  through  the  rotor  sideslip  angle 
to  align  controls  orthogonal  to  the  local  wind.  For¬ 
ward  rotor  cyclic  inputs  are  corrected  for  delta  three 
(63)  hinging  effects. 

7.  Coning  (aQ) ,  longitudinal  (ai) ,  and  lateral  (bi)  first- 
harmonic  flapping  coefficients  are  computed,  and  first- 
order  dynamic  lags  are  applied  to  account  for  rotor 
system  flapping  response. 

8.  Rotor  forces  and  moments  are  computed  using  either: 

a.  The  classical  equation  method,  with  forces  resolved 
perpendicular  and  parallel  to  the  shaft  normal 
plane  and  aligned  with  the  S.N.P.  wind  axis,  or 

b.  The  rotor  map  approach  with  thrust,  power  and  rotor 
propulsive  force  derived  from  tabular  isolated 
rotor  or  wind  tunnel  data  and  expressed  in  the 
rotor  wind  axis. 

9.  Fuselage  aerodynamic  forces  and  moments  are  calculated 
from  tabular  wind  tunnel  derived  data. 

10.  Rotor  aerodynamic  forces  and  moments  are  resolved 
through  rotor  sideslip  and  shaft  incidence  angles 
back  into  the  aircraft  body  axis.  Summation  of  these 
resolved  forces  and  moments  with  body  axis  fuselage 
forces  and  moments  is  performed,  and  to  these  values 
are  added  the  moments  imposed  upon  the  airframe  by 
individual  rotor  forces.  The  force  and  moment  sums 
are  inserted  into  the  equations  of  motion  to  again 
compute  the  airframe  acceleration  set. 

11.  Engine  governor  response  to  aerodynamic  rotor  power 
requirements  is  computed,  along  with  individual  rotor 
speed  degrees  of  freedom.  Rotor  shaft  windup  spring 
rates  are  accounted  for  in  the  governor  model. 

On  the  left-hand  side  of  Figure  5  are  shown  the  cockpit  and 
automatic  system  control  inputs.  The  automatic  systems  are 
not  covered  in  this  report  in  order  to  keep  the  rather  com¬ 
plicated  airframe  calculation  discussion  as  simple  as  possible 
for  clarity.  Reference  5  presents  a  comprehensive  explanation 
of  all  control  system  functions  along  with  the  appropriate 
modeling  equations.. 

Preceding  page  blank 
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TANDEM-ROTOR  MECHANICAL  CONTROL  SYSTEM  MODEL 


The  hybrid  mechanization  of  the  mechanical  and  automatic  con¬ 
trol  system  model  is  described  in  this  section.  It  includes 
the  347  conventional  mechanical  control  system. 

Figure  6  shows  a  functional  block  diagram  explaining  the  oper¬ 
ation  of  a  typical  tandem-rotor  helicopter  mechanical  control 
system.  Shown  on  the  left-hand  side  of  the  diagram  are  the 
lateral  (6S) ,  longitudinal  (6g)  ,  directional  (6r) ,  and  collec¬ 
tive  (Sc)  cockpit  control  inputs  made  by  the  pilot  or  gener¬ 
ated  analytically. 

Before  controls  are  mixed  mechanically,  the  high-authority 
differential  airspeed  hold  (DASH)  input  is  summed  with  the 
longitudinal  stick,  and  the  lower  directional  SAS  is  added  to 
the  pilot  pedal  input  as  shown  on  the  chart.  Control  mixing 
then  takes  place,  with  collective  pitch  inputs  (resulting  from 
63  and  6c  inputs)  going  to  each  head,  along  with  lateral 
cyclic  pitch  which  is  related  to  6r  and  6s  inputs.  Note  that 
the  collective  (9Cf#  9Cr)  and  differential  collective  (0Bf  and 
9Br)  controls  are  referenced  to  the  full-down  collective  posi¬ 
tion  on  each  rotor  head  (0Tp  and  0TR) . 

After  the  mechanical  control  mixing  is  accomplished,  SAS 
inputs  are  summed,  and  the  controls  are  passed  through  a 
second-order  upper  boost  model.  The  boost  actuators  are  des¬ 
ignated  as  being  either  swiveling  or  pivoting,  depending  upon 
their  primary  function  on  the  actual  aircraft  (note  the  S  or  P 
subscripts) .  For  both  rotors,  the  boost  actuators  operate 
either  together,  to  produce  collective  pitch,  or  differen¬ 
tially,  to  produce  lateral  cyclic  pitch. 

After  being  processed  through  the  upper  boost  system,  control 
inputs  are  "unmixed",  so  they  can  be  used  in  the  "classical" 
equations  or  with  the  rotor  maps.  Before  the  lateral  (A^c) 
controls  are  transformed  into  the  S.N.P.  wind  axis,  longitu¬ 
dinal  (B^c)  inputs  are  added.  Longitudinal  cyclic  pitch  for 
each  rotor  is  put  in  according  to  "q"  sensed  schedules  which 
vary  with  calibrated  airspeed.  Longitudinal  cyclic  may  be  put 
in  by  the  pilot  in  conjunction  with  the  AFCS. 

Longitudinal  cyclic  actuator  dynamics  as  shown  in  Figure  6  are 
modeled  with  a  simple  first-order  lag.  Because  of  the  rela¬ 
tively  high  frequencies  involved  in  the  SAS  actuators,  no 
dynamics  associated  with  their  operation  are  included. 

All  equations  describing  the  conversion  of  cockpit  controls  to 
the  equivalent  blade  pitch  angle,  control  mixing  equations  to 
accommodate  the  SAS  motion  reference,  and  the  final  conversion 
to  the  blade  pitch  angles  used  for  rotor  map  readings  are  de¬ 
tailed  in  Reference  5. 
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SMALL-MOTION  EXTERNAL  LOAD  MODEL 


The  aircraft  equations  of  motion  shown  schematically  in  Figure 
5  include  a  simplified  external  load  model  based  on  a  two-poir 
suspension  system.  The  modified  equations  are  shown  in  this 
section  together  with  load  aerodynamic  properties. 

The  simplified  external  load  model  is  based  on  three  additiona 
degrees  of  freedom  which  permit  longitudinal,  lateral  and  yaw 
motions  of  the  load,  all  measured  with  respect  to  the  airframe 
The  three  load  degrees  of  freedom  are  defined  by  cable  angles 
(see  Figure  7)  ySL,  XSL  and  vsl* 

The  three  additional  equations  of  motion  for  the  sling  load 
compute  angular  acceleration  of  the  cables,  i.e.,  ysL>  X‘sl» 
and  vsl-  Along  with  these  changes  are  terms  added  into  the 
original  six  airframe  equations  of  motion^ ref lectinq  coupling, 
which  is  expressed  as  a  function  of  ysL'  USL»  ^SL'  *SL»  VSL' 
and  vsl* 

The  external  load  equations  are  based  on  derivations  in  Refer¬ 
ence  5.  As  such,  they  are  the  simplest  ones  and  are  valid  for 
relatively  small  angular  motions  (15°  or  less)  of  the  two-poin 
parallel  cable  suspension.  A  general  external  load  model, 
based  on  total  force  approach,  which  includes  the  tether  re¬ 
straints  with  elastic  or  inextensible  cable,  is  described  in 
detail  in  Appendix  I. 

FULL-FORCE  EXTERNAL  LOAD  MODEL 


While  the  small-motion  external  load  analytical  model  describee 
in  Reference  5  is  adequate  for  most  hover  and  forward  flight 
simulation  tasks  in  which  cabl~  tensions  remain  nearly  con¬ 
stant,  it  is  difficult  to  adapt  to  special  transition  condi¬ 
tions  such  as  may  be  encountered  during  external  load  lift-off, 
from  both  cable  unloaded  through  partial  cable  tension,  to  the 
fully  airborne  state.  For  this  reason,  a  different  set  of 
equations  has  been  derived  to  accommodate  the  simulation  of 
these  special  flight  conditions.  These  equations  treat  the 
external  load  as  an  independent  rigid  body,  with  its  own  six 
degrees  of  freedom,  and  specify  a  set  of  helicopter  and  ex¬ 
ternal  load  forces  which  couple  the  two  bodies  by  means  of  two 
arbitrarily  located  elastic  or  inextensible  cables.  The  de¬ 
tailed  derivation,  based  on  the  Lagrangian  approach  applicable 
to  the  system  with  constraint  equations,  is  presented  in 
Appendix  I. 


32 


LONGITUDINAL  GEOMETRY 


X 


Figure  7.  Simplified  External  Load 
Model  Geometry. 
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Figure  8  illustrates  the  model  for  the  combined  external  load 
and  helicopter  system.  The  formulation  is  sufficiently  general 
to  accommodate  the  following  features: 

•  Arbitrary  location  of  two  hooks  on  the  aircraft  and 
attachment  points  on  the  load. 

•  Elastic  or  inextensible  cables. 

•  All  inertia  properties  of  the  external  load. 

In  the  preliminary  evaluation  of  the  elastic  cables  versus  the 
inextensible  cable  formulation,  it  was  established  that  the 
simulation  of  elastic  cables  will  be  simpler  for  two  basic 
reasons:  first,  in  the  case  of  inextensible  cables,  the  solu¬ 
tion  of  constraint  equations,  arising  from  forcing  the  cables 
to  be  infinitely  stiff  in  axial  direction,  will  require  more 
computation  to  restrain  the  motion,  including  the  computation 
of  linear  position,  velocities,  and  accelerations  along  the 
direction  of  cables;  second,  when  the  cables  become  transiently 
unstressed  (slack) ,  the  return  to  the  stressed  state  with  in¬ 
extensible  cables  will  require  the  computation  of  impulsive 
forces  which  must  be  applied  to  effect  a  stepwise  change  in 
linear  and  angular  (aircraft  and  load)  velocities  resulting 
from  sudden  tautness  of  the  cables.  The  elastic  cable  formu¬ 
lation  requires  only  a  continuous  tracking  of  cable  stretch 
magnitudes  and  rates.  The  transient  cable  slacking  and  stress¬ 
ing  will  be  automatically  accommodated  by  cable  elastic  forces 
alone. 
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Figure  8.  Full  Force  External  Load 
Model  (Side  View) . 


ANALYTICAL  STUDY  OF  ADVANCED  TANDEM- ROTOR 
HELICOPTER  (ATH)  HOVER  PERFORMANCE 


The  procedure  that  was  followed  in  analyzing  comparable  pre¬ 
cision  hover  performance  of  a  hover  hold  system  and  a  pilot 
was : 

1.  Investigate  unaugmented  helicopter  characteristics. 

2.  Improve  helicopter  handling  qualities  through  stability 
augmentation  loops. 

3.  Synthesize  hover  hold  system  and  estimate  hover  hold 
performance  in  turbulence. 

4.  With  hover  hold  system  off,  estimate  pilot-in-the-loop 
hover  hold  performance  and  pilot  workload  in  turbu¬ 
lence. 

UNAUGMENTED  AIRCRAFT 


To  represent  the  motion  of  the  rigid-body  helicopter,  a  set  of 
small-perturbation  differential  equations  was  used.  The  equa¬ 
tions,  as  presented  in  Figure  9,  represent  the  unaugmented 
helicopter,  separated  into  longitudinal-vertical  motion  and 
lateral-directional  motion.  This  separation  of  axes  was  main¬ 
tained  throughout  the  analysis.  In  later  sections,  these 
equations  were  modified  and  added  to,  to  include  the  effects 
of  gusts,  augmentation,  and  pilot  model.  The  equations  are  in 
first-order  form  suitable  for  state  matrix  evaluations. 

Long.'  tudinal-Vertical  Motion 

Figure  10  shows  the  effect  that  helicopter  size  has  on  the 
longitudinal-vertical  eigenvalues  (roots  of  the  characteristic 
equation) .  As  shown,  the  vertical,  or  plunging,  mode  (Zw)  re¬ 
mains  relatively  constant  with  increasing  size.  The  short¬ 
term  pitching  mode  (Mg)  also  remains  relatively  constant. 
However,  the  unstable  second-order  response,  which  is  reflec¬ 
tive  of  the  long-term  speed  and  attitude  changes,  shows  a 
slight  decrease  in  instability  and  a  slight  reduction  in 
frequency. 

Lateral-Directional  Motion 


Figure  11  shows  the  effect  that  helicopter  size  has  on  the 
lateral-directional  eigenvalues.  The  stable  aperiodic  mode 
that  corresponds  to  Nr  shows  little  effect  in  changing  size. 
The  other  aperiodic  mode  is  predominantly  roll  rate  damping 
(Lp)  and  shows  a  slight  decrease  in  time  to  half  amplitude 
with  increasing  size.  The  eigenvalue  at  the  origin  corres¬ 
ponds  to  the  neutral  stability  associated  with  heading 


36 


(Xq-Wo)q-g0 


M  |H 
+ 


+ 

> 


M 

tSJ 


in 


Z 

O 

H 

1 

O 

w 


z 

o 

H 

O 


H 

a 


w 

S 


• 

c 

o 

•H 

■P 


m 

0 

u> 

C 

0 

•H 

■P 

m 

3 

o1 

w 

■p 

M-t 

m 

P 

u 

p 

■H 

< 

P 

nJ 

<u 

c 

•H 

PJ 


ON 

0) 

3 

O' 

•rl 

IP 


37 


fa 

h 

fa 

O 

O 

w 

w 

< 

a 

w 

& 

38 


REAL  -  SEC  x 

Figure  10.  Longitudinal  System  Effect  of  Gross  Weight. 
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Figure  11.  Lateral-Directional  System  Effect  of  Gross  Weight 


changes.  The  slightly  unstable  second-order  pair  of  eigen¬ 
values  is  reflective  of  long-term  changes  in  lateral  velocity 
and  bank  angle.  A  slight  decrease  in  frequency  is  evident 
with  increasing  size. 

STABILITY  AUGMENTATION 


Stability  augmentation  in  the  form  of  feedback  loops  was  added 
to  the  control  system  to  improve  helicopter  responses.  To 
shorten  the  time  required  to  synthesize  a  good  system,  guide¬ 
lines  for  choosing  shaping  and  loop  gains  are  taken,  for  the 
most  part,  from  work  done  on  the  heavy  lift  helicopter  and  the 
Model  347. 

Loop  Closures 

To  show  the  effect  of  these  feedback  loops,  one  loop  at  a  time 
will  be  closed  with  various  loop  gains.  For  this  purpose  only 
the  50,000-pound  normal-gross-weight  configuration  will  be 
analyzed. 

Longitudinal-Vertical  System 

Figure  12  is  a  block  diagram  of  the  feedback  loops  added  to 
the  longitudinal-vertical  system  to  improve  pitch  rate  damping 
and  to  stabilize  the  long-term  second-order  response.  An  al¬ 
titude  rate  feedback  loop  into  collective  could  have  been 
added  at  this  time  to  improve  vertical  damping.  However,  the 
vertical  mode  is  almost  completely  uncoupled  from  the  longi¬ 
tudinal  modes,  and  thus  will  be  improved  when  the  precision 
hover  loops  are  added.  Closing  the  pitch  rate  feedback  loop 
first.  Figure  13  shows  the  effect  that  loop  gain  has  on  the 
locations  of  the  eigenvalues.  Increasing  loop  gain  causes  the 
eigenvalues  corresponding  to  the  pitch  rate  damping  mode  and 
feedback  lag  to  combine  and  form  a  second-order  coupled  re¬ 
sponse.  The  eigenvalues  representing  the  vertical  damping 
mode  and  feedback  washout  are  only  slightly  affected.  The 
long-term  unstable  pair  of  eigenvalues  moves  slightly  toward 
the  imaginary  axis,  but  remains  unstable  for  the  range  of 
gains  chosen.  The  maximum  gain  shown  is  the  final  value 
chosen  for  closure  of  this  loop,  and  all  subsequent  loop 
closures . 

With  the  rate  loop  closed.  Figure  14  shows  the  effect  of  clos¬ 
ing  the  pitch  attitude  loop.  As  loop  gain  increases,  the 
higher  frequency  pitch  mode  tends  to  decrease  in  frequency. 
However,  the  more  predominant  effect  of  this  loop  closure  is 
to  stabilize  the  long-term  response.  As  loop  gain  increases, 
the  unstable  pair  of  eigenvalues  crosses  the  imaginary  axis, 
approaches  the  real  axis,  and  splits  into  two  aperiodic  modes. 
As  gain  is  further  increased,  one  of  the  aperiodic  modes  de¬ 
creases  in  time  constant  while  the  other  increases.  The 
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Figure  12.  Longitudinal-Vertical  Stability  Augmentation. 


Figure  13.  Longitudinal  System  -  Effect  of  Pitch  Rate  Feedback 
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Figure  14.  Longitudinal  System  -  Effect  of  Pitch  Attitude  Feedback 


eigenvalue  corresponding  to  the  uncoupled  washout  time  con¬ 
stant  slightly  increases  its  time  constant.  Again,  as  with 
the  pitch  rate  feedback  loop,  pitch  attitude  feedback  into 
differential  collective  has  little  effect  on  the  vertical 
damping  mode. 

Lateral-Directional  System 


Figure  15  is  a  block  diagram  showing  the  selected  lateral- 
direccional  stability  augmentation  loops.  Figure  16  shows 
the  effect  of  closing  the  roll  rate  loop.  The  effect  is 
similar  to  that  of  the  pitch  rate  loop  closure,  in  that  the 
roll  rate  damping  (Lp)  mode  combines  with  the  feedback  lag  to 
form  a  second-order  coupled  response.  The  effect  of  roll  rate 
feedback  on  the  directional  modes  is  small.  The  effect  of 
this  loop  on  the  long-term  second-order  lateral  mode  is  a 
slight  improvement  in  the  stability. 

Figure  17  is  the  result  of  closing  the  roll  attitude  loop  with 
the  roll  rate  loop  closed.  As  the  gain  increases,  the  higher 
frequency  roll  rate  mode  decreases  in  damping.  The  long-term 
unstable  oscillation  becomes  stable,  approaches  the  real  axis, 
and  splits  into  two  aperiodic  modes.  As  loop  gain  is  further 
increased,  one  eigenvalue  becomes  more  negative,  reflecting 
the  tight  bank  angle  hold,  while  the  other  eigenvalue  combines 
with  the  washout  time  constant  to  form  a  lightly  damped  long¬ 
term,  low-frequency  oscillation.  The  directional  modes  are 
only  slightly  affected  by  this  roll  attitude  loop. 

Figure  18  is  the  result  of  adding  yaw  rate  feedback  to  the 
previously  closed  loops.  As  yaw  rate  feedback  gain  increases, 
only  the  mode  corresponding  to  Nr  is  affected.  The  time  to 
half  amplitude  of  this  mode  decreases  by  a  factor  of  about 
twenty. 

Figure  19  shows  the  effect  of  now  closing  the  yaw  attitude 
loop.  As  can  be  seen,  only  the  yaw  rate  mode  and  the  neutral¬ 
ly  stable  heading  mode  are  affected.  As  gain  increases,  these 
two  modes  combine  to  form  a  second-order  response. 

Figure  20  is  the  result  of  adding  an  inertial  lateral  velocity 
feedback  loop.  The  value  of  this  loop  can  be  seen  in  the  im¬ 
proved  long-term  pair  of  eigenvalues.  As  gain  increases,  this 
mode  increases  in  damping  ratio  until  it  splits  into  two  ape¬ 
riodic  modes. 

Effect  of  Helicopter  Size  With  SAS  ON 

With  all  the  loops  closed  as  shown  in  Figures  12  and  15,  the 
effect  of  increasing  the  helicopter  size  will  now  be  investi¬ 
gated. 
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Figure  15.  Lateral-Directional  Stability  Augmentation. 
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REAL  -  1/SEC 

Figure  16.  Lateral-Directional  System  Root  Locus  Effect  of  Roll  Rate  Feedback 
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Figure  17.  Lateral-Directional  System  Root  Locus  Effect  of  Roll 
Attitude  Feedback. 
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Longitudinal-Vertical  Motion 


Figure  21  is  a  plot  of  the  longitudinal-vertical  eigenvalues 
with  SAS  ON  for  four  different  gross  weights.  As  size  in¬ 
creases,  it  appears  that  the  higher  frequency  pitch  mode  ap¬ 
proaches  the  real  axis  and  splits  into  two  aperiodic  modes. 

As  gross  weight  continues  to  increase,  one  of  the  aperiodic 
eigenvalues  becomes  more  negative  while  the  other  moves  in  the 
direction  of  the  origin  until  it  combines  with  an  aperiodic 
pitch  attitude  mode  and  forms  a  second-order  response  at  a 
lower  frequency. 

Lateral-Directional  Motion 


Figure  22  shows  the  effect  of  increasing  helicopter  size  on 
the  lateral-directional  eigenvalues.  Both  of  the  second-order 
responses,  lateral  and  directional,  maintain  almost  constant 
damping,  but  experience  a  reduction  in  frequency. 

UNPILOTED  HOVER  HOLD  PERFORMANCE 

In  the  previous  paragraphs,  the  helicopter  was  augmented  by 
feedback  loops  to  improve  its  handling  qualities.  Even  though 
the  helicopter,  with  these  loops  closed,  is  now  stable,  there 
is  no  guarantee  that  it  will  return  to  its  original  position 
in  space  when  disturbed.  To  provide  position  hold,  position 
errors  can  be  used  to  feed  inputs  into  the  control  system. 

These  inputs  can  be  provided  either  automatically  through 
feedback  loops  or  by  the  pilot  through  optically  sensed  errors. 
This  section  deals  with  the  automatic  concept. 

Loop  Closures 

Feedback  loops  were  added  to  the  longitudinal-vertical  system 
and  the  lateral-directional  system  to  provide  automatic  hover 
hold  capability.  As  was  done  before,  the  effect  of  closing 
these  loops  one  at  a  time  was  analyzed  using  the  50,000-pound 
configuration. 

Longitudinal-Vertical  System 

Figure  23  is  a  block  diagram  of  the  selected  feedback  loops  to 
provide  longitudinal  and  vertical  position  hold.  Note  that 
the  longitudinal  signals  are  fed  back  into  longitudinal  cyclic 
and  not  differential  collective,  as  were  pitch  rate  and  pitch 
attitude.  This  is  done  because  longitudinal  cyclic  is  a  more 
effective  longitudinal  control,  with  less  pitch  coupling.  Al¬ 
though  the  longitudinal  and  lateral  positions  were  the  primary 
interest  in  this  report,  vertical  hold  capability  was  included 
for  completeness.  Since  it  has  been  shown  that  the  vertical 
motion  is  almost  completely  uncoupled  from  the  longitudinal 
motion,  any  improvement  in  vertical  position  hold  capability 
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Figure  21.  Longitudinal  System  -  Effect  of  Gross  Weight 
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Figure  22.  Lateral-Directional  System  -  Eigenvalues  Effect  of 
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Figure  23.  Longitudinal-Vertical  Precision  Hover  System. 


can  be  done  with  a  simple  one-degree-of-freedom  model,  inde¬ 
pendent  of  the  results  obtained  here  for  the  longitudinal  and 
lateral  system. 

Closing  the  inertial  longitudinal  velocity  loop  first.  Figure 
24  shows  the  result  of  various  feedback  gains.  As  can  be 
seen,  the  high-frequency  pitch  mode  is  affected  by  this  loop 
closure,  indicating  that  this  mode  is  now  strongly  coupled 
with  longitudinal  velocity.  Increasing  loop  gain  tends  to 
decrease  the  damping  in  this  mode  while  increasing  the  fre¬ 
quency.  This  loop  closure  also  causes  the  eigenvalue  associ¬ 
ated  with  velocity  damping  to  increase  negatively  until  it 
combines  with  the  pitch  attitude  mode  to  form  a  second-order 
response.  Further  increases  in  loop  gain  cause  this  mode  to 
increase  in  frequency  at  fairly  constant  damping  ratio.  This 
implies  that  the  rise  time  will  decrease  and  the  settling  time 
will  decrease  while  the  overshoot  remains  constant.  This  is 
a  desirable  condition  for  a  tight  and  rapid  control  system, 
provided  the  damping  level  is  sufficient  to  keep  overshoot 
within  tolerances.  The  eigenvalues  that  appear  at  the  origin 
represent  the  neutral  stability  of  altitude  and  longitudinal 
position.  These  loops  are  yet  to  be  closed. 

Now,  with  the  longitudinal  velocity  loop  closed,  the  longitu¬ 
dinal  position  loop  is  closed,  Figure  25.  Increasing  the  loop 
gain  moves  the  neutrally  stable  eigenvalue  from  the  origin  and 
drives  it  along  the  negative  real  axis,  thus  decreasing  its 
response  time  constant.  A  final  time  constant  of  about  1  sec¬ 
ond  is  reached.  One  eigenvalue,  representing  altitude,  still 
remains  at  the  origin. 

Figure  26  shows  the  effect  of  closing  the  altitude  rate  loop. 
The  effect  of  this  loop  closure  with  ever-increasing  gain  is 
to  drive  the  vertical  damping  time  constant  down,  with  the 
altitude  eigenvalue  remaining  at  the  origin. 

The  final  longitudinal-vertical  system  loop  closure  is  the 
altitude  loop  (Figure  27) .  The  effect  is  to  cause  the  aperi¬ 
odic  modes  representing  vertical  velocity  and  vertical  position 
to  combine  into  a  second-order  response.  Further  increase  in 
gain  causes  the  coupled  response  to  decrease  in  damping  and 
increase  in  frequency. 

Lateral  System 

Figure  28  is  a  block  diagram  of  the  selected  feedback  loops  to 
provide  lateral  position  hold.  Figure  29  shows  the  effect  of 
adding  a  lateral  velocity  feedback  through  a  washout.  The 
high-frequency  roll  rate  mode  increases  in  damping  while 
slightly  decreasing  in  damped  natural  frequency.  As  gain  in¬ 
creases,  two  lateral  aperiodic  modes  combine  to  form  a  second- 
order  lateral  response.  Further  gain  increases  cause  the 
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Figure  24.  Longitudinal  System  -  Effect  of  Longitudinal  Ground  Speed  Feedback. 
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Figure  26.  Longitudinal  System  -  Effect  of  Altitude  Rate  Feedback 


Figure  27.  Longitudinal  System  -  Effect  of  Altitude  Feedback 
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Figure  28.  Lateral  Precision  Hover  System. 
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Figure  29.  Lateral- Directional  System  -  Root  Locus  Effect  of 
Lateral  Velocity  Through  Washout  Feedback. 


damping  of  this  mode  to  decrease.  One  other  lateral  aperiodic 
mode  is  driven  toward  the  origin.  The  second-order  directional 
mode  is  unaffected. 

Figure  30  is  the  result  of  closing  the  lateral  position  loop 
with  all  other  loops  closed.  The  higher  frequency  roll  mode 
is  only  slightly  affected.  The  lower  frequency  roll-velocity 
mode  maintains  approximately  constant  damping  while  decreasing 
in  frequency.  The  neutrally  stable  eigenvalue  at  the  origin 
representing  lateral  position  is  pushed  out  the  negative  real 
axis,  consequently  quickening  its  response.  The  effective 
washout  time  constant  is  increased,  and  the  eigenvalue  is 
driven  closer  toward  the  origin.  The  second-order  directional 
mode  is  unaffected. 

Effect  of  Helicopter  Size  With  PHS  ON 

With  all  the  feedback  loops  representing  the  stability  augmen¬ 
tation  and  the  precision  hover  systems  closed,  the  effect  of 
helicopter  size  was  investigated. 

Longitudinal-Vertical  Motion 

Figure  31  shows  the  effect  of  increasing  helicopter  size  on 
the  longitudinal-vertical  responses.  Increasing  size  has 
almost  no  effect  on  the  second-order  vertical  response,  due 
to  the  little  change  in  Zw.  The  higher  frequency  second-order 
longitudinal  response  did  not  change  much  either.  This  points 
out  the  strong  dependence  of  this  mode  on  longitudinal  cyclic 
control  sensitivity,  which  is  relatively  invariant  with  size. 
However,  the  lower  frequency  second-order  longitudinal  re¬ 
sponse  shows  a  marked  decrease  in  frequency  at  almost  constant 
damping.  This  effect  is  due  to  the  dependence  of  this  mode  on 
DCP  control  sensitivity,  which  decreases  with  increasing  size, 
and  also  to  a  lesser  degree,  to  a  decrease  in  longitudinal 
cyclic  pitch  control  sensitivity. 

Lateral-Directional  Motion 


Figure  32  shows  the  effect  on  lateral-directional  responses 
due  to  increased  helicopter  size.  The  high-frequency  roll 
mode  increases  in  damping  and  decreases  in  frequency.  This 
effect  is  probably  due  to  increasing  Lp  and  decreasing  L6S 
with  increasing  size.  The  lower  frequency  second-order  roll- 
velocity  mode  shows  only  a  slight  decrease  in  damping.  This 
is  probably  due  to  the  strong  dependence  of  this  mode  on 
lateral  control  sensitivity  (Y6S),  which  remains  relatively 
constant.  The  second-order  directional  mode  shows  a  slight 
decrease  in  damping  and  a  definite  decrease  in  frequency, 
which  correspond  to  the  slight  decrease  in  Nr  and  definite 
decrease  in  NgR  with  size. 
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Figure  32.  Lateral-Directional  System  Eigenvalues 
Effect  of  Gross  Weight. 


Low- Speed  Gust  Model 

One  way  to  investigate  the  hold  performance  of  a  helicopter 
and  control  system  is  to  calculate  the  position  root-mean- 
square  (rms)  error  of  the  helicopter  subjected  to  random 
velocity  disturbances  (turbulence).  For  this  purpose,  a  one¬ 
dimensional  gust  model  was  chosen.  Random  U-gusts  were  used 
to  evaluate  the  longitudinal-vertical  motion,  while  random 
V-gusts  were  used  for  the  lateral-directional  motion. 

The  one-dimensional  gust  spectrum  chosen  for  this  analysis  is 

*  M  = 

“  +  “g 


where 


ajg  =  1 
9  2  Lt 

The  quantity  Lt  is  generally  considered  to  be  proportional  to 
altitude  at  low  altitudes  and  equal  to  1000  feet  at  altitudes 
above  1000  feet.  The  parameter  wg,  gust  break  frequency,  is 
chosen  to  be  tt/10,  agreeing  with  References  2  and  3.  The  value 
of  o g  is  chosen  to  be  5  ft/sec  for  most  of  the  study.  However, 
since  the  rms  of  any  response  is  proportional  to  crg,  the  result 
for  any  other  value  of  ag  is  easily  obtained.  Figure  33  is  a 
plot  of  this  gust  spectrum  and  shows  a  very  sharp  drop-off 
with  frequency. 

Hover  Hold  Performance  in  Turbulence 


Using  the  gust  model  just  described  and  the  fully  augmented 
helicopter,  the  performance  of  the  helicopter  in  turbulence 
was  investigated. 

Longitudinal-Vertical  Performance 


Figure  34  is  the  fully  augmented  (SAS  ON  and  PHS  ON)  longitu¬ 
dinal-vertical  system.  Figure  35  shows  the  modulus  squared  of 
the  transfer  function,  which  represents  the  Longitudinal  posi¬ 
tion  error  of  the  50,000-pound  helicopter  to  longitudinal  gust 
as  a  function  of  frequency.  This  plot  shows  the  maximum  value 
occurring  at  about  0.1  rad/sec.  Figure  36  is  the  longitudinal 
position  spectrum,  which  is  obtained  by  multiplying  the  modulus 
squared  of  Figure  35  and  the  gust  spectrum  of  Figure  33.  The 
predominant  effect  of  the  gust  spectrum  can  be  seen  in  the 
rapid  attenuation  with  frequency  of  Figure  36.  If  the  curve 
of  Figure  36  is  calculated  for  frequencies  between  minus  in¬ 
finity  and  plus  infinity,  then  integrated  between  these  limits, 
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Figure  33.  Gust  (Turbulence)  Power  Spectrum. 
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Figure  34.  Longitudinal-Vertical  System  -  Complete 
Block  Diagram. 
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Figure  35.  Modulus  of  Ratio  of  Longitudinal  Position  Error 
to  Longitudinal  Gust. 


and  finally  divided  by  2ir,  the  result  would  be  the  root-mean- 
square  response  of  position  error  due  to  the  gust  spectrum. 
However/  since  the  position  power  spectrum  is  greatly  attenu¬ 
ated  with  frequency/  an  adequate  approximation  to  the  rms  can 
be  obtained  by  truncating  the  integration  at  some  cutoff  fre¬ 
quency.  Figure  37  shows  the  result  of  varying  cutoff  frequency 
on  rms  position  error.  As  can  be  seen,  an  asymptotic  value  is 
rapidly  approached  above  1.0  rad/sec. 

A  measure  of  the  relative  comfort  and  discomfort  of  the  pilot 
is  the  rms  normal  acceleration  levels.  Figure  38  is  similar 
to  Figure  35  except  that  the  transfer  function  is  now  vertical 
acceleration  at  pilot  station  over  longitudinal  gust.  As  can 
be  seen,  the  maximum  value  is  obtained  at  a  much  higher  fre¬ 
quency  than  position  error.  This  is  due  to  the  pitch  acceler¬ 
ation  at  the  predominant  high-frequency  pitch  mode.  But  as 
shown  in  Figure  39,  the  power  spectrum  of  this  function  due  to 
the  gust  spectrum  is  quite  attenuated  above  4  rad/sec.  Figure 
40  shows  that  the  rms  pilot  acceleration  is  relatively  constant 
above  a  cutoff  frequency  of  18  rad/sec. 

Table  X  presents  the  rms  values  of  the  longitudinal  position 
error,  vertical  position  error,  pitch  attitude  error,  longitu¬ 
dinal  acceleration  at  c.g.,  vertical  acceleration  at  c.g. ,  and 
vertical  acceleration  at  pilot  station  as  a  function  of  heli¬ 
copter  size.  The  effect  of  size  is  minimal  on  all  parameters, 
with  the  greatest  effect  appearing  in  pilot  station  accelera¬ 
tion,  which  shows  a  slight  increase  in  size,  probably  due  to 
increasing  distance  between  c.g.  and  pilot  station. 

Lateral-Directional  Performance 


Figure  41  is  the  fully  augmented,  SAS  ON  and  PHS  ON,  lateral- 
directional  system.  Figure  42  is  the  lateral  position  error 
transfer  function  due  to  lateral  gust  for  the  50,000-pound 
configuration.  The  peak  value  for  this  function  occurs  at  a 
slightly  higher  frequency  than  did  the  longitudinal  position 
error,  and  does  not  drop  off  as  rapidly  with  frequency.  Fig¬ 
ure  43  is  the  corresponding  lateral  position  power  spectrum 
as  a  result  of  the  gust  spectrum  and  shows  a  slightly  shallow 
er  attenuating  slope  than  did  the  longitudinal  position  spec¬ 
trum.  As  a  result  of  this  shallower  slope,  the  rms  lateral 
position  error  approaches  a  higher  value  than  did  the  longitu 
dinal  position  rms  value,  Figure  44. 

Figure  45  is  the  modulus  squared  of  lateral  acceleration  at 
the  pilot  station.  The  first  peak  in  this  function  roughly 
occurs  at  the  lower  frequency  lateral  velocity  mode,  while 
the  second  and  smaller  peak  is  influenced  by  the  higher  fre¬ 
quency  lateral  mode.  Figure  46  is  the  result  of  imposing  the 
gust  spectrum  on  the  acceleration  modulus.  The  gust  spectrum 
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Figure  39.  Power  Spectrum  of  Vertical  Acceleration  at  Pilot 
Station  Due  to  Longitudinal  Gust. 
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Figure  40.  RMS  Vertical  Acceleration  at  Pilot  Station  Due  to 

Longitudinal  Gust  Versus  Integration  Cutoff  Frequency. 


TABLE  X. 

LONGITUDINAL-VERTICAL  MOTION  EFFECT  OF 

GROSS  WEIGHT  ON  RMS  RESPONSES  TO  GUST 

SAS  ON,  PHS  ON 

Gross  Weight 
(Lb) 

50,000 

180,000 

120,000 

200,000 

Xe  (Ft) 

0.03 

0.03 

0.03 

0.03 

Ze  (Ft) 

0.05 

0.05 

0.05 

0.05 

0e  (De9) 

0.13 

0.13 

0.13 

0.14 

ax 

c.g. 

(Ft/SecT) 

0.015 

0.015 

0.014 

0.014 

az 

c.g.  - 
(Ft/Sec^) 

0.025 

0.025 

0.025 

0.025 

a  z_ 

P  2 

(Ft/Secz) 

0.206 

0.207 

0.190 

0.215 

GUST  MODEL 

One-Dimensional  Longitudina 

*  -  y  W 

ug  (t)2  +  u)g2 

oug  =5.0  Ft/Sec 

tog  =  tt/10  Rad/Sec 

1  Spectrum 
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Figure  41.  Lateral-Directional  System  (Complete 
Block  Diagram) . 
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Figure  43.  Power  Spectrum  of  Lateral  Position  Erroi 
Due  to  Lateral  Gust. 
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^^■^ure  44.  RMS  Lateral  Position  Error  Due  to  Lateral  Gusts 
Versus  Integration  Cutoff  Frequency. 
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Figure  45.  Modulus  of  Ratio  of  Lateral  Acceleration  at  Pilot 
Station  to  Lateral  Gust. 
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Figure  46.  Power  Spectrum  of  Lateral  Acceleration  at  Pilot 
Station  Due  to  Lateral  Gusts. 


almost  eliminates  the  second  peak.  Finally,  the  effect  of 
cutoff  frequency  is  shown  in  Figure  47. 

Table  XI  presents  the  rms  values  for  the  lateral  position 
error,  bank  angle  error,  heading  error,  lateral  acceleration 
at  c.g.,  and  lateral  acceleration  at  pilot  station  as  a  func¬ 
tion  of  helicopter  size.  As  evident  from  the  table,  helicopter 
size  has  little  effect. 

PILOTED  HOVER  HOLD  PERFORMANCE 


The  ability  of  the  pilot  to  perform  the  precision  hover  task 
with  the  automatic  hold  system  off  is  investigated  in  this 
section.  A  linear  pilot  model  is  for  this  purpose. 

Pilot  Model 


The  pilot  transfer  function  model  conventionally  employed  in 
handling  quality  analyses  consists  of  the  general  form 

Kp  (TtS  +  1)  e~TS 
YP  =  (TnS  +  1)  (TjS  +  1) 

together  with  adjustment  rules  for  the  parameters  Kp,  Tl,  and 
Tj.  The  fixed  parameters r  and  T^  represent  reaction-time  and 
neuromuscular  lag.  General  rules  for  the  selection  of  these 
five  pilot  parameters  are  presented  in  Figure  48,  which  is 
reproduced  from  Reference  4. 

From  comments  made  in  Reference  5  and  preliminary  results  here, 
it  was  determined  that  the  adjustable  pilot  lag,  Tj ,  was  not 
desirable  for  position  loop  closures  in  the  crossover  fre¬ 
quency  range  of  interest.  Figure  49  is  representative  of  the 
effect  that  pilot  lag  has  on  closed-loop  stability.  Thus, 

Tj  =  0  is  used  during  the  analysis. 

The  actual  block  diagram  of  the  pilot  model  used  in  this  study 
is  shown  in  Figure  50.  The  transport  lag,  e-TS,  is  represent¬ 
ed  by  the  first-order  Pads ' s  approximation.  Values  for  the 
fixed  parameters  rand  tn  are  chosen  as  0.15  second  and  0.10 
second,  respectively. 

Control  Options  and  Pilot-Assist  Functions 

In  evaluating  the  pilot's  ability  to  perform  the  precision 
hover  task,  several  control  options  and  assist  functions  will 
be  considered.  Figure  51  is  a  block  diagram  showing  the  con¬ 
trol  options  and  assist  functions  in  the  longitudinal  axis. 
Three  longitudinal  options  will  be  studied.  Referring  to 
Figure  51, 
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TABLE  XI.  LATERAL-DIRECTIONAL  MOTION  EFFECT  OF 

GROSS  WEIGHT  ON  RMS  RESPONSES  TO  GUST 

SAS  ON,  PHS  ON 

Gross  Weight 

(Lb) 

50,000 

80,000 

120,000 

200,000 

Ye  (Ft) 

0.097 

0.096 

0.096 

0.095 

$e  (Deg) 

0.120 

0.120 

0.120 

0.120 

'I'e  (Deg) 

0.0115 

0.0115 

0.0115 

0.0172 

°YC.G.  0 

(Ft/SecT) 

0.054 

0.053 

0.053 

0.052 

%  2 
(Ft/SeO 

0.051 

0.050 

0.050 

0.048 

GUST  MODEL 

One-Dimensional  Lateral  Spectrum 

<Dv  =  2fa)g<7vg2 
g  w2+o>g2 

°vg  =  5.0  Ft/Sec 

wg  =  /10  Rad/Sec 
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Figure  49.  Effect  of  Adjustable  Pilot  Lag  on  Longitudinal 
Position  Loop  Closure  Stability. 


87 


(T  S+l) (-2S+1) 
P  (Tn  S+l)  (gS+1) 


Figure  50.  Pilot  Model  (Block  Diagram  Used  in  the  Analysis)  . 
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1.  Longitudinal  Control  Type  1:  longitudinal 

cyclic  control  only  (K6BIC  =  0.1, 

t<$BIC  =  °‘1'  Kob  =  0) 

2.  Longitudinal  Control  Type  2:  differential 

pitch  control  only  =  0, 

K6b  =  T(Sb  =  0) 

3.  Longitudinal  Control  Type  3:  longitudinal 

cyclic  plus  lagged  differential 
pitch  control  (K<sBIC  =  0.1, 

t6bic  =  ®  ^ '  ^<$b  =  ^ '  ^63  =  ^ 

Figure  52  shows  the  control  options  and  assist  functions  avail 
able  in  the  lateral  axis. 

1.  Lateral  Control  Type  1:  lateral  cyclic  control 

with  no  pickoffs  (K  Lcp  =  0) 

2.  Lateral  Control  Type  2:  lateral  cyclic  control 

with  roll  quickening  through  stick 
pickoff  (K  Lcp  =  1,  T^  =  0.8,  T^  =  3). 

In  addition  to  the  above  lateral  options,  one  further  possi¬ 
bility  will  be  investigated  in  an  effort  to  evaluate  the  use¬ 
fulness  of  velocity  feedback  in  aiding  the  pilot  to  perform 
his  task.  This  option  is  similar  to  option  1  except  that 
there  will  be  no  velocity  stabilization  through  the  SAS. 

3.  Lateral  Control  Type  0:  lateral  cyclic  control 

with  no  pickoffs  (KLCP  -  0)  and  no  velocity 

feedback  through  SAS  (G^  =  0) 


Pilot  Closure  of  the  Position  Loop 

In  this  study,  it  is  assumed  that  the  pilot  is  concerned  with 
one  axis  at  a  time,  either  longitudinal  or  lateral.  Thus  only 
single-axis  pilot  closures  will  be  addressed.  According  to 
Reference  4,  the  pilot  will  adjust  his  transfer  function  to 
perform  the  task  and  provide  good  closed-loop  stability.  The 
question  arises  as  to  which  takes  priority.  Reference  3  sug¬ 
gests  that  the  pilot  tends  to  minimize  pilot  rating  as  a 
weighted  function  of  performance  and  pilot  leads.  It  also 
suggests,  however,  that  a  minimum  20%  gain  margin  for  stability 
be  satisfied.  Reference  4  quotes  numbers  like  40  to  80  degrees 
phase  margin  and  closed- loop  damping  ratios  greater  than  0.350. 
Preliminary  results  showed  that  position  error  in  turbulence 
decreased  with  increasing  pilot  gain,  but  that  above  a  certain 
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Figure  52.  Pilot  Closure  of  Lateral  Position  Loop. 


value  the  stability  of  the  system  decreased.  Therefore,  it  is 
felt  that  the  pilot,  in  order  to  perform  the  task,  will  sacri¬ 
fice  as  much  stability  as  he  can.  Thus,  a  minimum  level  of 
0.350  for  closed-loop  damping  of  the  position  mode  will  be 
assumed  to  form  a  boundary  for  possible  combinations  of  pilot 
gain  and  lead  time  constants.  For  a  second-order  system,  this 
damping  level  corresponds  roughly  to  35-40  degrees  of  phase 
margin.  In  situations  where  it  is  felt  that  a  second-order 
response  is  not  a  good  approximation  for  the  position  response 
phase  margins  will  be  used  as  a  lower  stability  bound,  along 
with  the  requirement  that  no  mode  will  have  a  damping  ratio 
less  than  35%.  Once  these  stability  limits  on  pilot  gain  and 
lead  are  determined,  then  a  probable  range  of  predicted  values 
will  be  chosen  based  on  hold  capability  in  turbulence.  Pilot 
work  load,  in  terms  of  pilot  lead,  will  then  be  estimated. 

Longitudinal  Position 

Each  of  the  three  longitudinal  control  types  mentioned  above 
was  studied  using  the  50,000-  and  200,000-pound  configurations 
There  are  five  eigenvalues  of  interest  when  the  pilot  closes 
the  longitudinal  position  loop.  These  are  the  cubic,  which 
represents  the  predominant  pitch  response,  and  the  quadratic, 
which  represents  the  longitudinal  position  and  velocity  re¬ 
sponse.  For  low  pilot  gains  and  even  high  pilot  gains  under 
certain  control  types,  these  responses  can  be  considered 
distinct.  However,  with  the  200,000-pound  configuration  and 
control  type  3,  the  responses  become  highly  coupled.  For  the 
cases  where  the  pitch  response  becomes  lightly  damped,  the 
single  tracking  pilot  loop  hypothesis  becomes  doubtful.  In 
other  words,  when  pitch  excursions  become  significant  to  the 
pilot,  then  the  pilot's  inputs  to  the  control  system  are  going 
to  be  a  sensed  combination  of  pitch  and  position.  However, 
this  possibility  will  not  be  studied  as  it  is  felt  to  be 
beyond  the  scope  of  vhis  analysis.  In  these  cases  the  single¬ 
axis  approach  will  be  used  with  reservations. 

Longitudinal  Control  Type  1  -  Longitudinal  Cyclic 
Control  Only 

For  the  50,000-pound  configuration,  the  pitch  cubic  of  inter¬ 
est  is  a  complex  pair  at  relatively  high  frequency  and  a  real 
eigenvalue  with  a  longer  time  to  half  amplitude.  Pilot  clo¬ 
sure  of  the  position  loop  wa th  the  range  of  gains  studied  has 
little  effect  on  the  complex  pair  of  pitch  eigenvalues  and 
tends  to  drive  the  real  eigenvalue  toward  the  origin. 

Figure  53  is  a  root  locus  plot  of  the  complex  pair  of  eigen¬ 
values  representing  the  position-velocity  response.  Since  the 
aperiodic  pitch  mode  remains  significantly  quicker  (time  to 
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Figure  53.  Longitudinal  System  -  Effect  of  Pilot  Gain 
on  Closed-Loop  Stability. 
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half  amplitude  much  less)  than  the  position  mode  when  the  im¬ 
posed  stability  limits  are  approached,  it  can  be  assumed  that 
the  pitch  response  is  uncoupled  insofar  as  distracting  the 
pilot  from  his  precision  hover  <-ask  is  concerned. 

The  damping  of  the  position-velocity  mode  is  plotted  in  Figure 
54  as  a  function  of  pilot  gain  and  lead  time  constant.  For 
constant  pilot  lead,  damping  tends  to  increase  and  then  de¬ 
crease  as  pilot  gain  increases.  As  lead  increases,  the  shape 
of  the  damping  ratio  versus  pilot  gain  curve  shows  a  sharper 
peak  and  greater  attenuation  following  the  peak.  This  points 
out  a  definite  desirability  on  the  pilot's  part  to  provide  a 
little  lead  as  necessary  to  perform  the  task.  This  is  true 
because  the  higher  the  lead,  the  harder  the  pilot  would  have 
to  work  to  maintain  acceptable  bounds  on  his  gain  and  the  more 
the  stability  of  the  system  would  be  changing.  This  point  is 
further  emphasized  when  pilot  gain  is  plotted  against  pilot 
lead  for  constant  damping  levels.  Figure  55  is  such  a  plot 
for  the  50,000-pound  configuration  with  control  type  1.  As 
can  be  seen,  a  certain  amount  of  pilot  lead  is  needed  to  reach 
the  proposed  stability  limit.  The  pilot  will  choose  to  oper¬ 
ate  on  the  flat  portion  of  these  curves  if  he  can  satisfy  the 
task  there.  The  ability  of  the  pilot  to  perform  the  task  on 
any  portion  of  these  curves  will  be  considered  later  when  the 
hold  capability  of  this  configuration  is  analyzed  in  turbu¬ 
lence. 

For  the  200,000-pound  configuration,  the  pitch  cubic  is  a  low- 
frequency  complex  pair  and  a  large  negative  real  eigenvalue. 
Unlike  the  50,000-pound  configuration,  where  only  the  real 
eigenvalue  showed  significant  movement  with  pilot  gain,  the 
entire  pitch  cubic  is  affected  significantly  with  pilot  gain. 
Figure  56  shows  the  effect  of  increasing  pilot  gain  at  con¬ 
stant  pilot  leads  on  the  pitch  cubic  and  position-velocity 
quadratic.  As  gain  increases,  one  of  the  complex  pair  of 
pitch  eigenvalues  approaches  the  origin,  while  the  other  ap¬ 
proaches  the  large  negative  eigenvalue  that  is  moving  toward 
the  origin;  and  finally,  as  pilot  gain  further  increases,  the 
two  larger  eigenvalues  combine  to  form  a  complex  pair.  This 
complex  pair  continues  to  decrease  in  damping  and  increase  in 
frequency.  The  longitudinal  velocity-position  eigenvalues 
move  toward  each  other,  with  increasing  pilot  gain,  combine 
to  form  a  complex  pair,  and  then  increase  in  frequency  and 
decrease  in  damping  until  they  are  driven  unstable.  Figure  57 
is  the  damping  ratio  of  this  mode  versus  pilot  gain.  The  same 
trends  as  existed  for  the  50,000-pound  configuration  are  shown 
in  this  plot  for  the  200,000-pound  configuration.  Figure  58 
shows  the  damping  contours  for  the  200,000-pound  configuration 
and  shows  that  the  minimum  acceptable  pilot  lead  is  only 
slightly  higher  than  that  for  the  50,000-pound  configuration, 
indicating  that  the  position  response  for  this  type  of  control 
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Figure  54.  Longitudinal  System  -  Damping  of  Longitudinal 
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Figure  55.  Longitudinal  System  -  Damping  Boundaries  for 
Pilot  Gain  vs  Pilot  Lead  Time  Constant. 
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Figure  56.  Longitudinal  Root  Locus. 
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Figure  57.  Longitudinal  System  -  Damping  of  Longitudinal 
Position  Mode. 
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Figure  58.  Longitudinal  System  -  Damping  Boundaries 
for  Pilot  Gain  and  Lead. 


is  rather  invariant  with  size.  However,  the  pilot  might  find 
the  strong  pitch  coupling  in  the  200,000-pound  configuration 
undesirable . 

Longitudinal  Control  Type  2  -  Differential 
Collective  Pitch  Control  Only 

As  with  control  type  1,  the  closure  of  the  position  loop  with 
this  type  of  control  has  little  effect  on  the  high-frequency 
complex  pair  of  pitch  roots  for  the  50,000-pound  configuration. 
However,  as  shown  in  Figure  59,  the  aperiodic  pitch  mode  in¬ 
creases  in  time  constant  until  the  time  to  half  amplitude  of 
this  mode  is  about  the  same  as  that  for  the  position  mode. 

This  is  not  surprising  with  this  type  of  control,  since  the 
pilot  is  trying  to  generate  position  changes  with  primarily 
pitch  changes.  The  crossover  frequency  for  the  position  mode 
using  this  type  of  control  is  much  less  than  that  for  control 
type  1.  Figure  60  is  a  plot  of  damping  versus  pilot  gain  for 
control  type  2.  The  trend  is  similar  to  that  for  the  other 
control  type  and  gross  weight.  Figure  61  is  a  plot  of  the 
stability  limits.  When  compared  with  Figure  55,  the  necessary 
minimum  pilot  lead  seems  to  be  about  1.0  second  higher  with 
DCP  than  with  LCP,  implying  a  greater  work  load.  Also,  the 
maximum  acceptable  gain  is  slightly  less  with  DCP. 

With  the  200,000-pound  configuration,  an  aperiodic  pitch  mode 
tends  toward  the  origin  with  pilot  closure  of  the  position 
loop.  However,  the  aperiodic  mode  for  the  200,000-pound  con¬ 
figuration  is  formed  from  the  splitting  of  the  low-frequency 
second-order  pitch  mode,  as  shown  in  Figure  62.  The  charac¬ 
teristics  of  the  second-order  position  response  are  similar 
to  the  50,000-pound  case,  but  the  crossover  frequencies  are 
less.  As  Figure  63  shows,  the  damping  ratio  is  less  for  the 
200,000-pound  configuration  for  a  given  gain  and  lead.  Con¬ 
sequently,  as  Figure  64  shows,  the  minimum  allowable  pilot 
gain  is  less. 

Longitudinal  Control  Type  3  -  Longitudinal  Cyclic 
With  Lagged  Differential  Collective  Pitch  Control 

Figure  65  is  the  root  locus  plot  for  the  50,000-pound  config¬ 
uration  with  this  type  of  longitudinal  control.  The  high- 
frequency  second-order  pitch  mode  is  not  significantly  affect¬ 
ed  and  consequently  is  not  plotted.  The  aperiodic  pitch  re¬ 
sponse  falls  between  LCP  only  (type  1)  and  DCP  only  (type  2) 
as  discussed  in  the  above  paragraphs.  This  aperiodic  response 
is  closer  in  time  to  half  amplitude  to  the  position  mode  than 
with  LCP  only,  but  farther  than  with  DCP  only.  The  crossover 
frequencies  for  the  second-order  position  response  are  also 
between  those  for  the  LCP  only  and  DCP  only  types  of  control. 
Figure  66  shows  the  damping  ratio  variation  of  the  position 
mode  with  pilot  gain  and  lead  time  constant.  The  trend  of 
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Figure  59.  Longitudinal  Root  Locus. 
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Figure  60.  Longitudinal  System  -  Position  Modal  Damping 
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Figure  61.  Longitudinal  System  -  Damping  Ratio  Boundaries 
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Figure  62.  Longitudinal  Root  Locus. 
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Figure  63.  Longitudinal  System  -  Damping  of  Longitudinal 
Position  Mode. 
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Figure  64.  Longitudinal  System  -  Damping  Boundaries. 
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Figure  65.  Longitudinal  Root  Locus. 
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these  curves  is  again  similar  to  the  other  control  types. 

Figure  67  shows  the  stability  boundaries  for  varying  pilot 
gain  and  lead.  The  minimum  acceptable  pilot  lead  is  close  to 
that  for  DCP  only,  and  the  maximum  allowable  pilot  gain  is 
somewhat  less.  This  is  a  little  surprising  in  that  one  might 
expect  the  results  to  lie  somewhere  between  that  for  LCP  only 
and  DCP  only.  Obviously,  the  DCP,  even  though  lagged,  still 
dictates  a  stronger  pilot  lead.  However,  as  will  be  seen 
later,  the  performance  in  terms  of  rms  position  error  in  tur¬ 
bulence  is  more  like  LCP  only,  indicating  that  the  shorter 
rise  times  (for  constant  overshoot)  associated  with  LCP 
control  allow  for  a  tighter  hold  capability. 

Figure  68  is  the  root  locus  for  the  significant  eigenvalues  of 
the  200,000-pound  configuration  with  control  type  3.  This  plot 
shows  the  strong  effect  that  pilot  closure  of  the  position  loop 
has  on  the  pitch  cubic.  For  the  low  pilot  gains,  the  locus 
starts  out  looking  like  that  for  the  LCP  only  control.  How¬ 
ever,  as  pilot  gain  increases,  the  mode  which  started  out  as  a 
pitch  mode  is  driven  unstable  and  not  the  position  mode.  Thus, 
as  the  pilot  tries  to  minimize  the  position,  he  becomes  aware 
of  the  pitch  response,  until  eventually  his  concern  is  in  con¬ 
trolling  both.  At  this  time,  our  single-loop  approach  becomes 
dubious.  Since  in  the  high  gain  area,  the  second-order  re¬ 
sponse  that  started  as  the  combined  position-velocity  mode  is 
no  longer  a  good  representation  of  the  position  response,  the 
damping  of  this  second-order  response  will  not  determine  the 
stability  limits.  It  was  felt  that  phase  margin  versus  pilot 
gain  and  lead  offered  a  better  feel  for  the  stability  bounds. 
Figure  69  shows  the  effect  of  pilot  gain  and  lead  on  position 
error  phase  margin.  The  attenuation  of  the  phase  margin  with 
gain  is  greater  for  the  higher  leads.  Figure  70  shows  the 
relationship  between  pilot  gain  and  phase  margin  for  constant 
phase  margin.  A  lower  limit  of  40  degrees  will  be  used  for 
the  stability  boundary  for  this  configuration  and  control 
type  1.  The  probable  reason  that  the  minimum  lead  is  less 
stringent  for  the  200,000-pound  configuration  than  for  the 
50,000-pound  configuration  is  that  DCP  control  sensitivity 
decreases  at  200,000  pounds,  while  LCP  sensitivity  remains 
relatively  constant. 

Lateral  Position 


Each  of  the  three  lateral  control  types  mentioned  previously 
will  be  studied  using  the  50,000-  and  200,000-pound  configura 
tions.  There  are  six  eigenvalues  of  interest  when  the  pilot 
closes  the  lateral  position  loop.  These  are  the  eigenvalues 
representing  the  lateral  motion  of  the  helicopter.  The  three 
eigenvalues  representing  the  directional  motion  are  affected 
by  this  pilot  closure. 
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Figure  67.  Longitudinal  System  -  Damping  Boundaries. 
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Figure  68.  Longitudinal  Root  Locus. 
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Figure  69.  Longitudinal  System  -  Open  Loop  Phase  Margin 
of  Output  Position  to  Position  Error. 
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Lateral  Control  Type  0  -  Lateral  Cyclic  Control 
With  No  Pickoffs  and  No  Lateral  Velocity  Feedback 
Through  SAS 

This  control  type  is  simply  the  pilot  controlling  lateral  posi¬ 
tion  with  lateral  stick  motion,  with  no  stick  pickoffs  and  no 
lateral  velocity  damping  through  the  stability  augmentation 
system.  Any  additional  lateral  velocity  damping  over  and 
above  the  helicopter's  inherent  damping  must  be  provided  by 
the  pilot  through  his  lead  compensation.  When  comparing  this 
type  of  control  with  type  1  below,  the  benefit  of  velocity 
feedback  in  the  SAS  for  relief  of  pilot  work  load  is  obvious. 

For  the  50,000-pound  configuration,  the  eigenvalues  of  interest 
are  the  lateral  complex  pair  near  the  origin,  as  shown  in  Fig¬ 
ure  28  for  Gy  =  0,  and  the  eigenvalue  which  starts  out  at  the 
origin  with  the  position  loop  open.  Figure  71  shows  the  re¬ 
sult  of  position  loop  closure  on  the  complex  pair  of  eigen¬ 
values.  Again,  as  with  other  control  types,  the  greater  the 
pilot  lead  compensation,  the  higher  the  crossover  frequency. 

The  eigenvalue  which  originates  at  the  origin  is  driven  out 
the  negative  real  axis  until  it  approaches  the  value  of  the 
roll  attitude-lateral  velocity  washout.  Figure  72  represents 
the  variation  of  damping  ratio  of  the  complex  pair  of  eigen¬ 
values  as  a  function  of  pilot  gain  and  lead.  The  curves  show 
the  peakiness  which  results  from  increasing  pilot  lead.  The 
result  is  that  higher  pilot  leads  make  possible  higher  damping 
ratios  but  over  smaller  ranges  of  pilot  gain.  Therefore,  as 
before,  the  pilot  will  seek  a  lead  value  that  will  satisfy  his 
minimum  stability  level  and  no  more.  Figure  73  shows  the 
damping  ratio  contours  for  this  configuration.  Due  to  the 
slower  stability  changes  associated  with  the  flatter  portion 
and  the  smaller  work  load,  the  pilot  would  tend  to  operate  on 
this  portion  of  the  curves.  A  minimum  pilot  lead  of  1.5 
seconds  is  necessary  for  35%  damping. 

For  the  200,000-pound  configuration,  the  trends  are  the  same, 
as  shown  in  Figures  74,  75,  and  76.  However,  the  damping 
levels  attainable  at  the  same  pilot  leads  are  slightly  less 
for  the  200,000-pound  configuration.  Thus,  the  minimum  pilot 
lead  for  minimum  stability  would  be  slightly  greater  with  the 
200,000-pound  configuration,  as  shown  in  Figure  76. 

Lateral  Control  Type  1  -  Lateral  Cyclic 
Control  With  No  Pickoffs 


This  control  type  is  the  same  as  just  discussed  except  that 
the  stability  augmentation  system  is  providing  lateral  damping 
through  lateral  velocity  feedback.  The  advantage  of  this 
automatic  damping  is  immediately  obvious.  As  shown  in  Figures 
77,  78,  and  79  for  the  50,000-pound  configuration,  the  added 
damping  provided  by  the  control  system  allows  the  pilot  to 
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Figure  71.  Lateral  Root  Locus. 
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Figure  73.  Lateral  System  -  Damping  Ratio  Boundaries. 
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Figure  74.  Lateral  Root  Locus. 
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Figure  75.  Lateral  System  -  Lateral  Position  Modal  Damping. 
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Figure  76.  Lateral  System  -  Damping  Ratio  Boundaries 
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Figure  77 
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Figure  78.  Lateral  System. 


122 


LATERAL  -  DIRECTIONAL  SYSTEM 
DAMPING  BOUNDARIES 


0  0.1  0.2 


0.3  0.4  0.5 


PILOT  GAIN  (Kp)  -  IN. /FT 


Figure  79.  Lateral-Directional  System  Damping  Boundaries. 
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close  the  position  loop  with  lower  lead  compensation.  And, 
indeed,  he  can  close  the  loop  with  no  lead  compensation  at  all 
for  low  pilot  gains.  Figure  79  shows  the  variation  of  the 
stability  boundaries  with  pilot  gain  and  lead,  and  depending 
on  task  requirements,  a  minimum  pilot  lead  of  zero  is  possible. 

The  results  for  the  200,000-pound  configuration  are  similar  to 
those  for  the  50,000-pound  configuration  except  that  the  max¬ 
imum  pilot  gain  possible  for  a  given  pilot  lead  for  a  desired 
minimum  damping  ratio  is  less  for  the  200,000-pound  configu¬ 
ration.  Figures  80,  81,  and  82  present  the  results  for  the 
200,000-pound  configuration  with  this  control  type. 

Lateral  Control  Type  2  -  Lateral  Cyclic  Control 
With  Roll  Quickening  Through  Stick  Pickoff 

This  control  type  is  the  same  as  just  discussed  except  that 
roll  quickening  has  been  added  through  a  lateral  stick  pickoff 
through  a  lag  and  washout.  The  results  are  not  favorable.  As 
shown  in  Figures  83,  84,  and  85  for  the  50,000-pound  configu¬ 
ration,  much  lower  maximum  pilot  gains  are  possible  at  a  given 
pilot  lead  for  a  desired  stability  level.  Figure  85  shows  the 
further  disadvantage  of  rapid  stability  changes  with  pilot 
gain  changes.  The  200,000-pound  configuration  results,  Figures 
86,  87,  and  88,  are  similar  with  even  lower  maximum  pilot  gains 
possible. 

Piloted  Hover  Hold  Performance  in  Turbulence 


In  an  effort  to  restrict  the  acceptable  ranges  of  pilot  pa¬ 
rameters,  the  previous  section  showed  the  effect  on  stability 
of  these  parameters.  This  section  will  attempt  to  further 
pinpoint  the  pilot's  selection  of  parameters  by  relating  them 
to  the  precision  hover  task  in  turbulence. 

Longitudinal  Hold  Capability 

Using  the  three  longitudinal  control  types  discussed  and  the 
50,000-  and  200,000-pound  configurations,  longitudinal  position 
error  due  to  random  U-gusts  will  be  investigated. 

Figures  89  and  90  are  the  results  for  the  50,000-pound  config¬ 
uration  w-'th  control  type  1.  Figure  89  shows  the  rms  position 
error  versus  pilot  gain  for  constant  stability  levels,  where 
the  pilot  lead  is  adjusted  at  each  gain  to  provide  desired 
stability.  Figure  90  shows  the  rms  position  error  versus  pilot 
lead  for  constant  stability  levels,  where  the  pilot  gain  is 
adjusted  at  each  lead  to  provide  desired  stability.  For  a 
damping  ratio  of  0.35,  an  optimum  pilot  lead  of  about  2  sec¬ 
onds  and  an  optimum  pilot  gain  of  about  0.12  in. /ft  will  result 
in  an  rms  position  error  of  about  0.35  feet.  If  this  error  is 
assumed  to  be  sinusoidal  for  the  time  being,  then  an  rms  of 
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Figure  81.  Lateral-Directional  System  -  Lateral  Position 
Modal  Damping. 
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Figure  82.  Lateral-Directional  System  Damping  Boundaries. 
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Figure  83.  Lateral  Root  Locus. 
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Figure  84.  Lateral  System. 
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Figure  85.  Lateral-Directional  System  Damping  Boundaries. 
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Figure  86.  Lateral-Directional  System. 
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Figure  87.  Lateral-Directional  System. 
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Figure  88.  Lateral-Directional  System  Damping  Boundaries. 
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Figure  89.  Piloted  Longitudinal  Hold  Capability. 
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Figure  90.  Piloted  Longitudinal  Hold  Capability. 


0.35  would  mean  a  peak  value  of  +  0.5  feet,  which  for  a  pilot 
gain  of  0.12  in. /ft  would  mean  pilot  inputs  of  +  0.06  inch. 
Whether  this  sensitivity  is  below  the  pilot's  capability  will 
be  one  item  to  be  investigated  on  the  simulator.  If  control 
sensitivities  need  to  be  changed,  then  all  affected  feedback 
gains  will  have  to  be  reevaluated.  However,  most  of  the 
analysis  will  be  valid,  since  pilot  gain  times  control  sensi¬ 
tivity  is  usually  a  constant  multiplying  factor.  Thus,  a 
change  in  control  sensitivity  would  result  in  a  proportional 
change  in  pilot  gain.  Table  XII  presents  the  optimum  longitu¬ 
dinal  pilot  parameters  and  associated  errors  for  the  other 
control  types  and  gross  weights. 

Lateral  Hold  Capability 

Using  the  three  lateral  control  types  discussed  and  the 
50,000-  and  200,000-pound  configurations,  lateral  position 
error  due  to  random  V-gusts  will  be  investigated.  Table  XIII 
presents  the  optimum  lateral  pilot  parameters  and  associated 
lateral  position  errors  for  the  various  control  types  and 
configurations . 

Pilot  Rating 

In  previous  studies,  efforts  to  relate  pilot  work  load  and 
pilot  rating  to  stick  activity  were  based  on  rms  levels  and 
were  unsuccessful.  However,  References  3  and  4  associate 
pilot  lead  to  pilot  rating,  in  that  increased  compensation 
results  in  decreased  ratings.  Figure  91  is  produced  from 
Reference  4  and  shows  the  incremental  pilot  ratings  associated 
with  lead  compensation.  Tables  XII  and  XIII  also  include  the 
incremental  change  in  pilot  ratings  for  the  predicted  pilot 
leads . 
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TABLE  XII.  ADJUSTABLE  LONGITUDINAL  PILOT  PARAMETERS 
PREDICTED  "OPTItlAL"  OPERATION 


Gross 

Weight 

(Lb) 

Control 

Type 

Pilot 

Gain 

KP 

(In. /Ft) 

Pilot 

Lead 

tl 

(Sec) 

Pilot 

Lag 

(Sec) 

Position 

Error 

ox 

(Rms  Ft) 

Pilot 

Rating 

Change 

50,000 

1 

0.12 

2.0 

0 

0.35 

1.5 

200,000 

1 

0.11 

2.0 

0 

0.40 

1.5 

50,000 

2 

0.11 

3.0 

0 

1.20 

2.5 

200,000 

2 

0.08 

3.0 

0 

1.60 

2.5 

50,000 

3 

0.08 

2.5 

0 

0.45 

2.0 

200,000 

3 

0.11 

2.0 

0 

0.35  [ 

1.5 

(1)  Based 

on  Kp  and 

Tl  Yielding  Min 

o  at  0 

X 

35  Damping 

Ratio . 
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Figure  91.  Approximate  Variation  of  Incremental  Pilot  Rating 
With  Pilot  Lead  Time  Constant. 
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TABLE  XIII.  ADJUSTABLE  LATERAL  PILOT  PARAMETERS 
PREDICTED  ’’OPTIMAL"  OPERATION  (D 


Gross 

Weight 

(Lb) 

Control 

Type 

Pilot 

Gain 

KP 

(In. /Ft) 

Pilot 

Lead 

tl 

(Sec) 

Pilot 

Lag 

(Sec) 

Position 

Error 

ay 

(Rms  Ft) 

—  ^ 

Pilot 

Rating 

Change 

50,000 

0 

0.20 

D 

0.40 

1.0 

200,000 

0 

0.20 

B 

,r  /*’’ 

0.40 

1.0 

50,000 

1 

0.30 

0 

0.25 

0 

200,000 

1 

0.30 

1.0 

0 

0.25 

0 

50,000 

2 

0.08 

0.5 

0 

0.60 

0 

200,000 

2 

0.08 

0.5 

0 

0.60 

0 

(1)  Based  on  and  TL  Yielding  Min  a ^  at  0.35  Damping  Ratio. 
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BASIC  SIMULATOR  MODEL  VALIDATION  APPROACH 


A  summary  of  the  Boeing-Vertol  flight  test  data  is  outlined  in 
this  section,  and  the  approach  to  the  hybrid  model  validation 
is  described.  The  full-flight  envelope  hybrid  simulation 
mathematical  model  used  in  the  Boeing-Vertol  simulation  facil¬ 
ity  was  validated  against  flight  test  data  taken  on  the  "wing¬ 
less"  347  helicopter.  The  hybrid  model  configuration  for 
these  validation  test  runs  included  the  "basic"  347  (airframe 
and  control  system)  and  the  same  vehicle  with  an  external 
sling  load  attached.  A  set  of  hybrid  computer  responses  was 
then  generated  for  the  future  use  in  the  validation  of  the 
Northrop  simulator  program.  The  complete  data  validation 
package,  which  is  fully  detailed  in  Reference  5,  was  submitted 
to  the  Northrop  Corporation  as  scheduled  in  the  contract. 

This  section  represents  only  a  short  synopsis  of  the  validation 
approach  and  Boeing-Vertol  simulator  validation  results. 

MATHEMATICAL  MODEL  VALIDATION  APPROACH 


The  basic  approach  to  obtaining  a  validated  simulator  model 
involved  the  simulation  on  the  Boeing-Vertol  hybrid  computer 
of  the  full-envelope  math  model  (Reference  5) ,  followed  by 
validation  of  this  model  against  piloted  Model  347  flight  test 
data  based  on  specified  pilot  inputs  or  electronically  gener¬ 
ated  SAS  inputs.  After  this  simulator  model  was  adequately 
validated,  a  set  of  response  data  was  generated  based  on  an 
analog  representation  of  the  347  control  system.  An  identical 
computational  setup  was  then  used  to  program  the  full-flight 
envelope  model  at  the  Northrop  facility,  and  the  total  new 
program  was  checked  out  against  the  validated  Boeing-Vertol 
simulation  results. 

BASIC  MODEL  347  AIRCRAFT  FLIGHT  TEST  AND  HYBRID 
SIMULATION  MODEL  COMPARISON 


AFCS  Configuration 

The  347  airframe  and  control  system  models  were  first  vali¬ 
dated  with  static  and  dynamic  response  data  generated  from 
theory.  Mechanization  of  the  automatic  flight  control  system 
functions  for  this  program  checkout  phase  was  accomplished 
with  an  analog  computer. 

The  analog  FCS  configuration  reflected  a  design  version  of  the 
system  that  was  later  optimized  through  testing  on  the  347 
flight  vehicle.  Changes  in  system  gains,  time  constants,  and 
variations  in  shaping  were  made  during  the  flight  test  opti¬ 
mization. 
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Prior  to  comparing  hybrid  response  data  with  flight  test  re¬ 
sults,  it  was  necessary  to  modify  the  347  control  system  model 
to  match  the  final  "optimized"  flight  test  installation.  This 
was  done  by  programming  the  AFCS  digitally  in  the  360/44  com¬ 
puter  (through  the  use  of  the  VECEX  modeling  tool  discussed  in 
Reference  5)  . 

The  control  system  model  was  synthesized  from  the  larger  347 
demonstrator  digital  AFCS  by  "enabling"  the  appropriate  rate, 
attitude,  pickoff,  and  synchronizer  loops  (with  the  proper  347 
gains,  time  constants,  etc.).  Details  of  this  digital  control 
system  mechanization  are  presented  in  Reference  5. 

Correlation  Results 


Flight  test  data  utilized  for  model  verification  included  a 
series  of  pulses  and  steps  put  in  by  the  pilot  (at  various 
airspeeds)  in  the  longitudinal,  lateral,  and  directional  axes. 
Selected  pulses  were  also  applied  through  the  aircraft  SAS. 

In  the  hybrid  lab,  these  test-derived  cockpit  control  deflec¬ 
tions  were  modeled  as  a  function  of  time  with  a  BCA  mini¬ 
computer.  BCA  outputs  were  trunked  to  the  appropriate  "control 
system"  analog  console  for  mixing  and  summation  with  model  SAS 
inputs. 

Figures  92,  93,  and  94  illustrate  typical  correlations  of 
flight  test  and  hybrid  results.  Figures  92  and  93  depict  an 
unaugmented  control  input  at  medium  cruise  speed  in  level 
flight.  This  plot  indicates  clearly  (through  the  good  corre¬ 
lation  of  rate  and  attitude  response  data)  that  the  airframe 
model  is  adequate  for  the  purposes  intenued.  The  SAS  ON  hover 
response  presented  in  Figure  94  also  shows  good  agreement  with 
expected  results.  Other  runs  made  for  correlation  with  flight 
data  are  discussed  in  detail  in  Reference  5. 
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DIRECTIONAL  RT.  PULSE  PILOT  INPUT 
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1 - 1 - 1 - 1 - 1 - 1 - ' - 1 - 1 - ' - r 

0  1  2  3  4  5  6 

TIME  -  SEC 


HYBRID  DATA _ 

DATE  AUG.  31,  1972 
RUN  49 

VECEX  347  AUTO  CONTROL 
SYSTEM  WITH  ROLL  RATE 
FEEDBACK  COMPENSATION 
INCLUDED 


TOGW  45,660  LB 
TOGG  9  IN.  AFT 
IAS  80  KT 
SAS  &  DASH  OFF 
RPM  220 
ALT(Hd)  2,800  FT 


FLT  TEST  DATA 

FLIGHT  428 
RUN  24 

AIRCRAFT  B-164 


FLIGHT  TEST  CASE  3 


Figure  92.  Hybrid  Model  Validation  With  347  Flight 
Test  Data. 
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LATERAL  RT.  PULSE  PILOT  INPUT 
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G  = . 5"/DEG 
9 

DATE  AUG.  31,  1972 
RUN  NO.  48 


T.O.  G.W.  45,660  LB 
T.O  G.G.  9  IN.  AFT 
IAS  85  KT 
SAS  &  DASH  OFF 
RPM  220 

ALT.  (Hd)  2,800  FT 


FLT.  TEST  DATA 

FLIGHT  420 
RUN  21 

AIRCRAFT  B-164 
FLT.  TEST  CASE  2 


Figure  93.  Hybrid  Model  Validation  With  347  Flight 
Test  Data. 
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Longitudinal  Pulse  Pitch  Down  Pilot  Input 
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CASE  5 


Figure  94.  Hybrid  Model  Validation  With  347  Flight 
Test  Data. 
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GROUND-BASED  SIMULATOR  EVALUATION  OF  HOVER  PERFORMANCE 


The  hover  performance  evaluation  was  accomplished  at  the 
Northrop  Corporation,  Aircraft  Division,  flight  simulation 
facility  utilizing  the  large  motion-base  simulator.  The  math¬ 
ematical  model  programmed  on  this  simulator  is  described  in 
FULL-FLIGHT  ENVELOPE  MATHEMATICAL  MODEL  and  in  Reference  1. 

The  precision  hover  system  (PHS)  loops  were  added  to  the  basic 
helicopter  model,  and  an  oscilloscope  display  of  the  helicop¬ 
ter  position  error  was  installed  in  the  cockpit  dashboard. 

All  piloted  evaluation  was  performed  by  a  Boeing-Vertol  pilot 
with  extensive  experience  in  experimental  and  production  heli¬ 
copters,  both  with  and  without  external  load. 

PRECISION  HOVER  SYSTEM  MECHANIZATION 

The  PHS  was  mechanized  on  an  analog  computer  with  the  differ¬ 
ential  position  error  computed  from  integration  of  lateral  and 
longitudinal  helicopter  velocities  which  were  resolved  into 
the  true  horizontal  plane.  These  integrated  velocities  were 
used  directly  as  the  signals  coming  from  a  synchronized  posi¬ 
tion  sensor.  Appropriate  logic  signals,  actuated  by  the  stick 
detent  and  PHS  switches,  were  added  for  position  synchroniza¬ 
tion  whenever  the  pilot  entered  the  control  loop.  A  continu¬ 
ous  velocity  integration  signal  drives  the  cockpit  display  of 
helicopter  position  with  respect  to  a  ground-fixed  reference. 

FEEDBACK  LOOPS  AND  SYNCHRONIZING  LOGIC 


The  PHS  feedback  loop  mechanization  is  shown  in  Figures  95, 

96,  and  97.  These  loops  are  the  additional  ones  to  those  de¬ 
signed  for  the  Model  347  simulation  (see  FULL-FLIGHT  ENVELOPE 
MATHEMATICAL  MODEL  and  Reference  5  for  details)  .  The  most 
significant  feature  of  the  longitudinal  PHS  is  in  the  longi¬ 
tudinal  cyclic  pitch  (LCP)  control  loops.  The  longitudinal 
velocity  and  position  errors  are  fed  back  directly  into  the 
LCP  control  on  both  rotors.  The  pilot  inputs  have  the  options 
of  direct  differential  collective  pitch  (DCP)  and/or  LCP  con¬ 
trol.  Both  of  these  options  were  used  in  the  piloted  simu¬ 
lator  evaluation.  The  ground  speed  and  the  ground  position 
error  feedback  loops  are  similar  for  all  three  (longitudinal, 
lateral,  and  vertical)  axes.  The  ground  position  reference 
synchronization  is  achieved  by  controlling  the  operate/reset 
modes  of  the  ground-speed  integrators  with  switching  logic  L8 
shown  in  Figure  97.  Lateral  ground  speed  and  position  error 
actuate  the  lateral  cyclic  pitch  control;  the  vertical  velocity 
and  the  altitude  error  move  the  collective  pitch  control.  The 
horizontal  position  loops  (lateral  and  longitudinal)  go  into  a 
synchronizing  mode  whenever  the  pilot  displaces  the  cyclic 
stick  out  of  detent.  The  altitude  hold  remains  closed  unless 
the  collective  stick  magnetic  brake  is  released. 
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I 

L8A - 1 

L8A  =  1:  PERATE  (INTEGRATE) 

LS\  =  0:  ESET 


Figure  95.  Longitudinal  Precision  Hover  System. 
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fiCAFCS 


L8C  =  1:  OPERATE 
L8C  =  0:  RESET 


Figure  96.  Vertical  Hold  System. 
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<S 

SAFCS 


L8B  =  1:  OPERATE 
L8B  =  0:  RESET 


Figure  97.  Lateral  Precision  Hover  System  (Sheet  1  of  2). 
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L3  PRECISION  HOVER  POSITION 
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L37 

■|  AFCS  ENGAGED  -  ALL  AXIS 


Figure  97.  Lateral  Precision  Hover  System  (Sheet  2  of  2). 
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COCKPIT  DISPLAY  FUNCTIONS 


A  scope  display  of  this  aircraft  position  with  respect  to  a 
ground-fixed  reference  was  provided  in  the  cockpit.  The  dis¬ 
play  mechanization  included  two  options:  an  "inside-out" 
view,  in  which  the  motion  of  the  scope  displayed  "T"  repre¬ 
sented  the  view  of  the  ground-fixed  reference  from  the  air¬ 
craft  (i.e.,  the  scope  is  fixed  to  the  aircraft) ;  and  an 
"outside-in"  view"  in  which  the  displayed  "T"  represented  the 
motion  of  the  aircraft  as  seen  by  an  observer  stationed  on  a 
ground-fixed  platform  above  the  aircraft  (i.e.,  the  scope  is 
fixed  to  the  ground) .  These  two  options  were  programmed  in 
deference  to  the  fact  that  not  all  pilots  like  the  same  view 
for  simulator  work.  Some  pilots  prefer  to  work  with  the 
"outside-in"  views,  and  thus  are  able  to  interpret  the  motions 
on  the  display  screen  as  the  actual  aircraft  behavior,  while 
other  pilots  find  the  opposite  arrangement  easier  to  follow. 
The  pilot  who  was  "flying"  the  simulator  in  this  evaluation 
used  the  "outside-in"  view  exclusively. 

The  basic  variables  which  were  used  to  generate  the  two  dis¬ 
play  functions  are  defined  in  Figure  98.  The  transformations 
for  the  two  options  are  defined  in  Figure  99. 

PILOT-ASSIST  MODES  AND  CONTROL  TYPES 


The  control  types  described  in  the  Piloted  Hover  Hold  Per¬ 
formance  section  were  modeled  on  the  simulator  for  a  compara¬ 
tive  evaluation  in  piloted  hover  performance.  In  addition, 
two  more  control  types  were  simulated  in  each  of  the  two  axes. 
In  the  longitudinal  axes  (see  Table  XIV) ,  these  are  control 
types  4  and  5;  in  the  lateral  axes  (see  Table  XV),  these  are 
control  types  3  and  4.  The  lateral  control  types  3  and  4  and 
the  longitudinal  control  type  4  were  added  in  order  to  provide 
the  pilot  with  ground  pattern  tracking  capability  over  large 
distances  (+  50  feet  or  more) .  The  longitudinal  control  type 
5  was  added  in  order  to  evaluate  the  effect  of  ground  speed 
damping  in  both  axes,  it  was  anticipated  that  this  type  of 
control,  obtainable  with  type  5  in  longitudinal  and  type  1  in 
lateral  axes,  should  be  one  of  the  best  available  to  the  pilots 
for  precise  tracking  tasks  over  small  distances  (+5.0  feet 
around  the  target) .  In  addition  to  the  control  types  shown  in 
Tables  XIV  and  XV,  the  pilot  could  maneuver  within  the  tight 
PHS  control  with  automatic  position  synchronizing  on  the 
cyclic  control  stick.  This  control  type  was  termed  the  PHS 
control.  As  such,  it  represented  the  baseline  control  type, 
and  all  Cooper-Harper  rating  changes  for  other  controls  were 
compared  with  respect  to  the  PHS  control  type. 
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NORTH  4> 


(y x)  =  COMPUTED  E-W  COORDINATE,  POSITIVE  EAST 
(l|>  )  =  HEADING  ANGLE,  POSITIVE  CLOCKWISE  FROM  NORTH 

(X ,Y)  =  COORDINATE  SYSTEM  FIXED  TO  GROUND  (INERTIAL) 
(x,y)  =  AIRCRAFT  BODY  COORDINATES 

Figure  98.  Scope  Display  Variable  Definition. 
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TABLE  XIV.  LONGITUDINAL  CONTROL  TYPES 
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X/4B  =  5.0  f ps/in. 


XSC'  ^SC'  ^sc  = 
SCOPE  DRIVING 


NORTH 


EAST 


Figure  99.  Cockpit  Scope  Display  Driving  Functions. 
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TABLE  XV.  LATERAL  CONTROL  TYPES 
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SIMULATION  PROCEDURES  AND  RESULTS 


The  simulator  evaluation  of  the  hover  performance  in  turbu¬ 
lence  was  done  only  for  the  50,000-pound  configuration.  The 
aircraft  performance  runs  without  the  pilot  control  were  made 
with  the  full  motion,  PHS-ON,  while  the  pilot  kept  his  hands 
off  the  controls.  The  piloted  performance  runs  were  made  for 
several  selected  control  types  at  nominal  turbulence  levels  of 
0.0,  5.0,  10.0,  and  15.0  feet  per  second  rms.  For  presenta¬ 
tion  purposes,  the  actual  turbulence  levels  were  used,  which 
were  calculated  from  the  digital  records  of  the  data  sampled 
at  five  samples  per  second,  in  addition,  a  tracking  task  (CTR) 
was  performed  by  the  pilot  with  each  of  the  selected  control 
types,  in  which  the  pilot  traced  a  closed  triangular  path  on 
the  ground  with  the  aircraft.  An  X-Y  plotter  was  used  to  show 
how  the  pilot  was  able  to  perform  these  tasks  at  various  tur¬ 
bulence  levels  with  different  control  types. 

A  typical  hover  performance  run  was  conducted  by  first  letting 
the  aircraft  attain  a  trim  configuration,  and  then  switching 
into  the  operate  mode  with  the  cockpit  motion.  Several  trial 
runs  were  then  obtained  (and  recorded  digitally)  to  check  on 
the  natural  drift  characteristics  of  the  cockpit  display  func¬ 
tions.  The  unpiloted  hover  performance  runs  then  followed,  in 
which  the  turbulence  level  was  varied  while  the  pilot  kept  his 
hands  off  the  controls.  The  pilot's  comments  with  respect  to 
the  acceleration  level  acceptability  and  general  realism  of 
the  cockpit  motion  were  noted.  Then  the  pilot  was  given  the 
hover  hold  task  (HH)  and  the  closed  path  tracking  task  (CTR) 
for  the  different  control  types  at  various  turbulence  levels. 
The  aircraft  response  data  was  recorded  at  the  selected  in¬ 
tervals,  and  the  pilot's  comments  regarding  Cooper-Harper 
rating  and  general  aircraft  response  were  noted  in  the  log. 

The  full  log  of  all  runs  is  shown  in  Appendix  III. 

Hover  Hold  Tasks 

The  piloted  hover  hold  tasks  are  described  in  this  section. 
These  tasks  were  performed  by  requiring  the  pilot  to  hold  the 
displayed  aircraft  position  (display  option  2,  Figure  99)  over 
the  fixed  ground  point.  This  particular  pilot  preferred  to 
work  with  "outside-in"  display  exclusively.  The  pilot  was 
given  the  hover  hold  task  with  the  baseline  PHS  control  capa¬ 
bility  (synchronized  position  hold)  and  with  control  type 
combinations  5/1  (long,  type  5/lateral  type  1),  1/2,  1/0,  2/0, 
and  1/1.  On  the  initial  runs,  the  pilot  found  control  combi¬ 
nations  1/0  and  2/0  hard  to  work  with,  apparently  because 
these  controls  did  not  hr.ve  either  lateral  or  longitudinal 
velocity  feedback.  With  the  additional  complication  of  small 
display  motions  (+5.0  feet  in  both  directions),  the  pilot  had 
difficulty  in  his"”attempts  to  remain  within  the  display  limi¬ 
tations  in  higher  than  5.0-fps  turbulence  levels.  The  set 
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of  X-Y  plots  shown  here  is  therefore  limited  to  control  types 
5/1,  2/1,  and  the  PHS  mode.  The  reduced  data  for  other  con¬ 
trol  combinations  is  discussed  in  the  Simulator  Evaluation  and 
Analytical  Results  section. 

The  baseline  PHS  control  type  hover  hold  capabilities  in  tur¬ 
bulence  are  shown  in  Figure  100.  This  control  was  considered 
to  be  the  best  by  the  pilot,  with  an  approximate  Cooper-Harper 
rating  of  3.0.  On  the  statistical  basis,  this  control  shows 
the  capability  to  hold  the  hover  position  with  about  4.0 
inches  (standard  duration)  in  15.0-fps  nominal  turbulence 
level,  in  practical  terms  this  means  that  in  high  turbulence 
the  pilot  will  be  able  to  hold  the  aircraft  within  4.0  inches 
of  the  fixed  ground  position  68.3%  of  the  total  duration  of 
the  task,  and  about  95.6%  of  the  task  duration  the  aircraft 
position  error  will  not  be  greater  than  8.0  inches  in  either 
direction.  Figures  101  and  102  compare  the  piloted  performance 
at  5.0  fps  nominal  turbulence  level  with  control  type  combi¬ 
nations  5/1  and  1/2  (longitudinal/lateral).  The  advantage  of 
the  longitudinal  velocity  feedback  (type  5/1)  is  clearly  visi¬ 
ble  in  the  significantly  better  position  hold  shown  in  Figure 
101  compared  with  Figure  102  (control  type  1/2) .  A  similar 
advantage  can  be  discerned  at  10.0  fps  nominal  turbulence 
levels  shown  in  Figures  103  and  104.  Performance  with  control 
combination  5/1  at  15.0  and  20.0  fps  turbulence  levels  is  dis¬ 
played  in  Figures  105  and  106.  Comparison  of  these  control 
combinations  with  the  baseline  PHS  control  mode  performance  of 
Figure  100  clearly  indicates  the  powerful  effect  of  a  good 
precision  hover  system  on  hover  hold  performance  in  turbulence. 

Hover  Tracking  Tasks 

In  addition  to  the  hover  hold  tasks  on  the  simulator,  the  pilot 
performed  closed-path  tracking  tasks  (CTR) ,  in  which  he  was 
instructed  to  fly  the  aircraft  around  a  triangular  path,  with 
4.0-ft  rectangular  sections.  This  section  describes  the 
pilot's  performance  of  the  CTR  tasks  at  different  turbulence 
levels  using  control  type  combinations  5/1,  1/2,  and  the 
baseline  PHS  control  mode. 

The  basic  tracking  capability  of  the  three  control  types  (in 
calm  air)  is  compared  in  Figures  107,  108,  and  109.  The  base¬ 
line  PHS  mode  (Figure  107)  emerges  as  the  best  control  type, 
control  type  5/1  (Figure  108)  rates  second,  and  control  type 
1/2  rates  third.  In  terms  of  Cooper-Harper  scale,  this  pilot 
rated  the  baseline  PHS  and  control  type  5/1  as  3.0.  control 
type  1/2  was  rated  as  5.0  at  the  full  design  longitudinal 
control  sensitivity,  and  4.0  when  the  longitudinal  stick  sen¬ 
sitivity  was  reduced  to  one-half  of  the  design  setting.  An 
improvement  in  longitudinal  tracking  capability  with  reduced 
stick  sensitivity  is  evident  when  the  traces  in  Figure  111 
are  compared  with  those  of  Figure  110.  The  tracking  capability 
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(b)  10.0  FT/SEC  NOMINAL  TURB.  LEVEL 


(c)  15.0  FT/SEC  NOMINAL  TURB.  LEVEL 


Figure  100.  Simulator  Results:  Piloted  Hover  Hold 
Performance  in  Turbulence. 


hrJM- 


3  JULY  1973 
RUN  NO.  14 


(WEST) 


5.0  FPS  NOMINAL  TURB. LEVEL 


COMPUTED  RESULTS: 

X  =  -0.158  FT 
OX  =  0.427  FT 
y  =  -0.130  FT 
Oy  =  0.515  FT 

U  =  -1.66  FPS 
ocL  =  4.45  FPS 
Vg  =  -2.56  FPS 
avg  =  8.56  FPS 


X  (SOUTH) 

SCALE:  1.0  FT/IN. 


Figure  101.  Piloted  Hover  Hold  Control  Type  5/1. 
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Figure  102.  Piloted  Hover  Hold  Control  Type  1/2. 


10  FPS  TURB. (NOMINAL  LEVEL) 


3  JULY  1973 
RUN  NO.  19 


Y  -* 
(WEST) 


SCALE:  1.0  FT/IN. 


X 

(SOUTH) 


COMPUTED  RESULTS 


X  =  -0.166 
—  0.683 

y  -  0.217 

Oy  =  0.645 


Ug  *  2.77  FPS 

aUg=  14.05  FPS 
Vg  =  3.59  FPS 

a„  =  13.95  FPS 

Vg 


Figure  103.  Piloted  Hover  Hold  Control  Type  5/1. 


NOMINAL  TURB. LEVEL:  10.0  FPS 
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L6 


3  JULY  1973 
RUN  NO.  5 


Figure  107. 


Piloted  CTR  Task  Control  Type  phs 


3  JULY  1973 
RUN  NO.  13 


SCALE:  1.0  FT/IN. 


Figure  108.  Piloted  CTR 


Task  Control  Type  5/1. 


3  JULY  1973 
RUN  NO.  7 


Figure  109.  Piloted  CTR  Task  Control  Type 


NO  TURBULENCE 


1/2  (Full  Gains). 


Figure  110.  Piloted  CTR  Task  Control 
Longitudinal  Gain) . 


3  JULY  1973 
RUN  NO.  9 

NO  TURBULENCE 
SCALE:  1.0  FT/IN. 


Type  1/2  (With  1/2  X 
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CTR  Task  Control  Type  PHS. 


is  further  compared  with  PHS  and  control  type  5/1  at  higher 
turbulence  levels  in  Figures  112  and  113. 

Effects  of  Limited  Simulator  Motion 


The  simulator  motion  cues  as  they  might  affect  the  piloted 
hover  hold  performance  were  explored  only  in  a  cursory  manner. 
Two  piloted  hover  hold  tasks  were  performed  with  and  without 
the  cockpit  motion.  Figure  114  shows  a  piloted  hover  hold 
run,  fixed  base  in  a  5.0-fps  turbulence  level;  Figure  115 
shows  the  same  run  repeated  with  the  moving  cockpit.  A  com¬ 
parison  in  absolute  performance  terms  clearly  indicates  that 
cockpit  motion  cues  helped  the  pilot  in  performing  the  same 
task.  Figure  116  shows  another  fixed-base  hover  hold  run  at 
10.0  fps  turbulence  level.  A  comparable  run  with  cockpit 
motion  is  shown  in  Figure  117,  but  this  run  was  performed  at 
a  later  date  when  the  pilot  evidently  had  learned  the  system. 
In  any  case,  judging  from  these  runs  and  the  pilot's  comments, 
good  cockpit  motion  significantly  improves  the  pilot's  ability 
to  perform  comparable  tasks. 
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Figure  115.  Piloted  HH  Task  Control  Type  5/1. 


CONTROL  TYPE:  5/1 
PILOT  TASK:  HH 
TURB.  LEVEL:  10.0  FPS 


FIXED-BASE  RUN 


27  JUNE  1973 
RUN  NO.  35 
(DIGITAL  FILE  28) 


Figure  116.  Piloted  Hover  Hold  Task  Control  Type  5/1. 
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Figure  117.  Piloted  Hover  Hold  Task 


3  JULY  1973 
RUN  NO.  15 

10  FPS  NOMINAL 
TURBULENCE  LEVEL 


COMPUTED  RESULTS 

X  =  -0.356  FT 
Ox  =  0.645  FT 

y  -  0.103  FT 

ay  =  0.339  FT 

Uq  =  1.96  FPS 

aSg  -  10.81  FPS 

Vq  -  -0.95  FPS 

aVg  =  6*87  FPS 


Control  Type  5/1. 
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GROUND-BASED  SIMULATOR  EVALUATION  AND  ANALYTICAL  RESULTS 

The  simulator  results  at  50,000  pounds  gross  weight  are  com¬ 
pared  in  this  section  with  the  analytical  predictions  described 
in  ANALYTICAL  STUDY  OF  ADVANCED  TANDEM-ROTOR  HELICOPTER  (ATH) 
HOVER  PERFORMANCE.  The  unpiloted  (SAS-ON,  PHS-ON)  simulator 
results  for  hover  hold  performance  are  compared  with  the  the¬ 
oretical  predictions  at  turbulence  levels  up  to  about  20.0 
fps.  The  comparison  is  made  in  terms  of  standard  duration 
position  errors,  with  the  inherent  assumption  that  over  a  long 
duration  of  the  hover  hold  task,  the  average  drift  is  either 
small  or  nonexistent.  The  piloted  results  for  hover  position 
hold  are  compared  on  the  same  basis,  with  the  additional  com¬ 
parison  of  how  the  actual  pilot  behaves  compared  to  the 
theoretical,  "optimal",  pilot  controller.  The  "optimal"  pilot 
controller  is  defined  here  as  that  controller  which  performs 
with  minimum  at  minimum  position  error  at  a  total  system 
damping  ratio  of  0.35.  As  will  be  seen  from  the  simulator 
data,  the  real  pilot  does  not  always  behave  in  this  "optimal" 
fashion,  but  tends  to  adjust  his  equivalent  transfer  function 
with  learning  to  the  point  of  sacrificing  stability  for  good 
performance,  or  searching  for  a  comfortable  work  load  at  the 
sacrifice  of  performance.  The  pilot's  Cooper-Harper  rating 
changes  are  also  compared  for  the  selected  control  type  com¬ 
binations. 

COMPARISON  OF  UNPILOTED  HOVER  PERFORMANCE 

The  aircraft  performances  with  all  automatic  systems  operational 
(SAS-ON,  PHS-ON),  but  without  active  pilot  control,  are  com¬ 
pared  directly  in  Figure  118.  The  solid  lines  in  this  figure 
represent  the  theoretical,  linear  system  estimates  for  longi¬ 
tudinal  and  lateral-directional  aircraft  position  errors  as  a 
function  of  rms  turbulence  levels.  Since  the  theoretical  sys¬ 
tems  are  linear,  the  errors  vary  linearly  with  the  turbulence 
levels.  The  simulator  data  points  are  plotted  directly  on 
this  figure.  However,  it  must  be  pointed  out  that  the  simu¬ 
lator  points  have  not  been  corrected  for  the  drift  due  to  a 
parallel  digital  integration  of  position  error  displayed,  com¬ 
pared  to  the  analog  feedback  error.  As  such,  the  simulator 
results  will  inherently  show  significantly  larger  errors,  in 
a  somewhat  random  fashion,  with  larger  durations  at  higher 
turbulence  levels.  As  seen  from  Figure  118,  the  simulator 
results  compare  well  with  the  theoretical  predictions  at 
normal  turbulence  levels  (5.0  fps),  but  diverge  somewhat  at 
high  turbulence  levels. 

COMPARISON  OF  PILOTED  HOvER  PERFORMANCE 

Ti.*s  theoretical  basis  for  piloted  performance  prediction  dur¬ 
ing  precise  hover  position  hold  (see  Piloted  Hover  Hold  Per¬ 
formance  section  for  details)  reduces  the  pilot  transfer 
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Figure  118.  Unpiloted  Hover  Hold  Performance  Analysis  and 
Simulator  Results. 


function  variation  to  two  basic  parameters:  the  pilot  antic¬ 
ipation  factor  (pilot  lead  time,  Tj,)  #  representing  his  work 
load,  and  the  pilot  gain  factor  (Kp) ,  denoting  the  pilot's 
stick  motion  per  unit  of  observed  position  error,  in  the 
analysis,  these  two  factors  determine  the  actual  closed-loop 
performance  of  the  pilot  and  the  aircraft.  The  theoretical 
pilot  is  treated  as  an  "optimal"  controller  in  the  sense  that 
the  pilot  will  try  to  adjust  his  lead  time  and  gain  factor  so 
as  to  minimize  the  aircraft  position  while  still  maintaining 
the  total  (closed-ioop  pilot  and  aircraft)  system  damping  in 
the  principal  control  mode  at  a  comfortable  level,  in  most  of 
the  references  dealing  with  the  study  of  pilot  dynamics,  this 
total  system  damping  level  is  found  to  be  in  the  neighborhood 
of  0.35  damping  ratio.  For  a  given  control  type,  the  analyti¬ 
cal  results  may  be  summarized  in  the  form  of  typical  maps  such 
as  the  one  shown  schematically  in  Figure  119.  This  map  de¬ 
fines  the  total  pilot-aircraft,  closod-loop  system  character¬ 
istics  plotted  in  terms  of  constant  damping  ratio  and  position 
error  (for  a  given  level  of  turbulence)  contours  plotted  versus 
pilot  lead  time  and  gain  factor.  The  solid  curves  indicate 
the  closed-loop  system  damping  ratio  contours;  the  dashed 
curves  denote  the  closed-loop  system  position  errors  for  a 
constant  turbulence  level.  For  this  particular  control  type, 
the  "optimal"  pilot  controller  should  operate  in  the  desig¬ 
nated  cross-hatched  area  where  the  position  error  is  about 
0.2  ft  rms  and  the  total  system  damping  ratio  is  0.35.  Oper¬ 
ation  in  the  upper  right-hand  side  of  this  area  denotes  regions 
of  good  performance  with  high  work  load  (high  pilot  lead  time)  . 
Operations  in  the  lower  right-hand  region  represent  good  per¬ 
formance  with  lower  than  the  "optimal"  damping  ratio.  On  the 
simulator,  the  real  pilot  will  adjust  his  transfer  function 
as  he  learns  the  system,  until  he  strikes  a  balance  between  a 
satisfactory  performance  and  the  total  system  damping.  In¬ 
tuitively,  one  might  expect  the  pilot  to  work  in  the  regions 
of  high  damping  when  he  begins  to  learn  the  system,  and  to 
sacrifice  the  damping  levels  as  he  becomes  more  confident  of 
his  control  capability.  Most  of  the  initial  runs  with  various 
ntrol  types  tend  to  indicate  the  real  pilot  operation  in  the 
ea  denoted  as  "observed  pilot  performance".  This  is  the 
region  in  which  the  total  system  damping  remains  constant  with 
pilot  gain  adjustment,  so  the  pilot  is  able  to  concentrate  on 
error  minimization  with  little  change  in  work  load. 

The  reduction  and  interpretation  of  piloted  simulator  hover 
hold  performance  follow  this  general  line  of  thought.  Once  a 
particular  form  of  pilot  transfer  function  is  assumed,  time 
histories  of  position  errors  and  stick  motion  can  be  used  to 
calculate  pilot  parameters.  By  assuming  that  the  pilot's 
neuromuscular  lag  and  reaction  time  delay  are  constant,  then 
pilot  gain  and  lead  can  be  calculated  using  a  least-squares 
error  technique.  With  time  histories  generated  from  tne 
piloted  simulation  effort  at  Northrop,  the  data  on  the 
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left-hand  side  of  Tables  XVI  and  XVII  was  calculated.  It 
shows  pilot  gain  and  lead,  root-mean-square  gust,  and  root- 
mean-square  position  error.  For  comparison  with  theory,  each 
combination  of  pilot  gain,  pilot  lead,  and  gust  level  was 
inputted  into  the  analytical  model,  and  root-mean-square  posi¬ 
tion  error  and  modal  damping  were  calculated.  Also,  the 
theoretically  predicted  optimums  for  pilot  gain  and  lead  for 
each  configuration  are  shown. 

Lateral  Piloted  Performance 


At  this  time,  certain  results  pertaining  to  the  lateral  axis 
have  been  reduced  (see  Table  XVI)  .  The  agreement  between 
simulator  position  errors  and  the  corresponding  analytically 
calculated  position  errors  is,  in  general,  quite  good,  pre¬ 
dicted  optimum  pilot  gains  also  show  good  agreement.  The 
optimum  predicted  leads  are  not  bad  for  lateral  control  types 
0  and  1;  however,  the  predicted  lead  for  control  type  2  is 
considerably  less  than  those  calculated  from  simulator  results. 
Certain  conclusions  will  be  drawn  from  the  data  presented  in 
Table  XVI,  but  it  must  be  remembered  that  not  all  the  simu¬ 
lator  data  has  been  reduced. 

Looking  at  the  one  data  point  for  the  0  type  lateral  control, 
the  gain  and  lead  show  good  agreement  with  predicted  optimums. 
However,  only  one  point  has  been  reduced,  and  certain  ques¬ 
tions  arise.  Referring  to  Figure  73,  the  analytical  results 
show  that  the  pilot  has  chosen  a  gain  below  the  optimum,  which 
yields  only  about  25%  modal  damping  with  the  lead  chosen. 
Speculating  on  one  data  point,  a  preliminary  conclusion  may 
be  that  the  pilot  was  still  feeling  the  system  out  and  had  not 
yet  determined  its  optimums.  Indeed,  the  pilot  did  not  fly 
this  configuration  for  very  long.  More  familiarization  with 
this  type  of  control  might  have  led  to  increased  performance, 
but  the  work  load  would  have  remained  relatively  high. 

The  data  points  for  the  lateral  control  type  1  show  a  gain 
variation  which  encompasses  the  predicted  optimum.  The  pilot 
leads  are,  in  general,  slightly  less  than  the  predicted  opti¬ 
mum.  Referring  to  Figure  120,  it  shows  that  stability  is 
relatively  constant  over  a  rather  wide  gain  range  for  damping 
ratios  of  35%  and  less.  Therefore,  as  the  pilot  becomes  more 
confident  with  the  control  type,  he  is  able  to  increase  his 
gain  and  decrease  the  stability  level  in  an  effort  to  minimize 
position  error. 

The  simulator-generated  position  errors  as  a  function  of  turbu¬ 
lence  level  for  control  type  1  are  compared  in  Figure  121.  it 
shows  that  the  assumption  of  pilot  performance  in  the  region  of 
constant  system  damping  at  a  damping  ratio  slightly  less  than 
0.35  gives  an  excellent  correlation  with  the  simulator  results, 
although  this  is  not  the  "optimal"  controller  operation  point. 
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TABLE  XVI.  LATERAL  PILOTED  HOVER  PERFORMANCE  COMPARISON  OF  SIMULATOR- 
MEASURED  AND  ANALYTICALLY  CALCULATED  DATA 
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TABLE  XVII.  LONGITUDINAL  PILOTED  HOVER  PERFORMANCE  COMPARI 
SIMULATOR-MEASURED  AND  ANALYTICALLY  CALCULATED 
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SAS  -  ON 
PHS  -  OFF 
50,000-LB  CONFIG. 
GUST  LEVEL  =5.0  FPS 


PILOTED  SIMULATION  POINTS 


Ye  =  LAT.  POSN.  ERROR,  FT  RMS 


PILOT  GAIN  (Kp)  -  IN./FT 

Figure  120.  Lateral-Directional  Piloted  Performance 
Analysis  and  Simulation  Results  Control 
Type  1. 
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Finally,  looking  at  control  type  2,  the  simulator  pilot  gains 
show  good  agreement  with  the  predicted  optimum.  However,  the 
pilot  leads  generated  by  the  pilot  on  the  simulator  are  con¬ 
siderably  higher  than  the  predicted  optimum.  Referring  to 
Figure  122,  the  stability  sensitivity  of  the  position  mode  is 
high  with  respect  to  pilot  gain,  in  other  words,  small  changes 
in  pilot  gain  cause  rather  large  stability  changes,  as  measured 
by  modal  damping,  in  predicting  an  optimum,  it  was  assumed 
that  such  a  system  would  encourage  the  pilot  to  sacrifice  some 
performance  in  order  to  minimize  lead;  thus,  instead  of  pre¬ 
dicting  a  pilot  gain  of  0.16  and  pilot  lead  of  2.25,  a  pilot 
gain  of  0.08  and  a  lead  of  0.5  were  chosen,  in  comparison 
with  the  simulator  results,  the  gain  predicted  looked  good, 
but  the  pilot  chose  to  create  a  larger  lead.  This  higher  lead 
at  this  gain  yielded  higher  damping  ratios,  45%  to  50%  criti¬ 
cal.  With  more  familiarization  with  the  control  type,  the 
pilot  may  have  relaxed  his  stability  criteria  and  either  in¬ 
creased  his  gain  for  increased  performance,  decreased  his  lead 
for  decreased  work  load,  or  both. 

Longitudinal  Piloted  Performance 

The  tabulated  results  for  longitudinal  piloted  performance  are 
shown  in  Table  XVII.  The  analytical  predictions  based  on  the 
"optimal"  controller  operation  show  the  same  trends  as  the 
simulator  results;  the  best  performance  is  obtainable  with 
control  type  5  and  the  poorest  with  control  type  2  (pure  DCP 
control)  ,  with  control  type  1  in  the  middle.  The  computed 
position  errors  based  on  actual  simulator  results  for  pilot 
lead  and  pilot  gain  are  generally  in  agreement. 

The  comparative  results  for  control  type  5  are  examined  in  de¬ 
tail  in  Figures  123  and  124.  Several  data  points  are  available 
for  this  control  type,  so  the  following  comparison  is  valid. 

The  position  errors  in  terms  of  turbulence  levels  (Figure  123) 
indicate  that  the  "optimal"  pilot  controller  predicts  the  best 
operation  capability,  which  is  much  better  than  the  average 
attained  by  the  pilot  from  the  simulator  evaluation.  In  Figure 
124,  this  discrepancy  is  revealed  to  be  the  result  of  the  dif¬ 
ferent  real-pilot  operation  from  the  theoretical  "optimal". 

On  the  average,  the  real  pilot  performed  in  the  area  of  good 
damping  but  with  very  little  work  load  (Tl  =  0.4  second)  .  it 
is  conceivable  that  the  pilot  chose  this  region  because  the 
total  system's  damping  level  gradient  is  relatively  steep  for 
this  control  type,  but  good  performance  can  be  obtained  even 
with  no  pilot  lead  requirement. 

Pilot  Rating  Summary 


The  pilot  ratings  for  different  control  types  obtained  from  the 
simulator  evaluation  are  based  on  only  one  pilot's  opinion,  who 
had  a  limited  time  available  for  evaluating  the  different 
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Figure  122.  Lateral-Directional  Piloted  Performance  Analysis  and 
Simulation  Results  Control  Type  2. 
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systems.  In  addition,  the  limited  scope  display  (+  5.0)  made 
the  pilot  apprehensive  during  operations  with  the  relatively 
poor  control  configurations,  because  of  the  tendency  for  the 
display  to  run  rapidly  off  the  scope.  However,  several  im¬ 
portant  points  may  be  noted  upon  examination  of  the  pilot's 
Cooper-Harper  ratings  shown  in  Table  XVIII.  First,  the  base¬ 
line  PHS  and  the  control  combination  5/1  were  rated  by  the 
pilot  as  equal.  This  is  in  agreement  with  the  analytical  pre¬ 
dictions.  Second,  the  least  acceptable  control  was  the  2/0 
combination,  representing  pure  DCP  control  combined  with  roll 
acceleration.  This  is  also  in  agreement  with  the  theory,  al¬ 
though  the  absolute  magnitude  of  predicted  rating  changes  is 
considerably  smaller  than  that  obtained  from  the  real  pilot. 

As  was  indicated  before,  the  main  apparent  reason  was  dis¬ 
agreement  in  the  absolute  magnitudes,  which  stems  from  the 
fact  that  the  pilot  does  not  behave  exactly  like  the  theoreti¬ 
cal  optimal  controller.  In  addition,  he  becomes  concerned 
when  the  display  has  limited  displacement  capabilities  and 
tends  to  rate  the  poorer  systems  lower,  sometimes  uncontrol¬ 
lable,  because  once  he  loses  the  display,  he  is  not  able  to 
locate  it  and  bring  it  back  within  the  visible  scope  area. 

GENERAL  ASSESSMENT  OF  GROUND-BASED  SIMULATOR  RESULTS 

Conclusions  Based  on  Present  Data 

Two  basic  conclusions  were  reached.  The  first  is  that  velocity 
feedback  is  desirable  in  that  it  relaxes  the  need  for  the  pilot 
to  create  large  leads.  The  velocity  feedback  thus  reduces  the 
pilot's  work  load  and  improves  his  ratings  of  the  system.  The 
other  conclusion  deals  with  the  part  that  system  stability 
plays  in  influencing  the  pilot's  choice  in  compensation.  This 
relates  also  to  the  pilot's  learning  curve,  if  the  stability 
of  the  system  changes  rapidly  with  pilot  gain  or  lead,  then 
the  pilot  will  tend  to  choose  his  parameters  such  as  to  close 
the  position  loop  with  rather  significant  damping  (data  re¬ 
duced  to  date  shows  damping  ratios  in  the  order  of  50%) .  On 
the  other  hand,  if  stability  remains  fairly  constant  with 
pilot  gain,  for  instance,  then  the  pilot  will  tend  to  operate 
at  a  lower  stability  level  if  it  yields  greater  performance. 
Also,  as  the  pilot  becomes  more  confident  in  the  system,  he 
will  tend  to  push  the  system  to  its  limits  if  performance  is 
improved.  As  stated  earlier,  the  task  is  of  primary  concern 
in  selection  of  pilot  compensation,  with  stability  character¬ 
istics  playing  a  key  secondary  role. 
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TABLE  XVIII.  PILOT  RATING  OF  CONTROL  TYPES 
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CONCLUSIONS 


The  following  conclusions  are  derived  from  the  studies  con¬ 
ducted  during  this  program: 

1.  The  tandem-rotor  helicopter  configurations  in  the 
VHLH  class  can  be  configured  so  that  their  hover 
hold  performance  is  similar  to  (or  directly  scalable 
from)  the  existing  MLH  class,  in  particular,  the 
control  system  configurations,  required  for  precise 
automatic  hover  position  hold  and  pilot  control, 
can  be  synthesized  from  the  MLH  applicable  designs 
and  related  experience.  The  singular  exception  to 
this  rule  might  be  in  the  rotor-induced  coupled 
elastic  fuselage  and  blade  lead-lag  modes.  Special 
filtering  might  be  required  to  preclude  an  unde¬ 
sirable  aggravation  of  these  low  damping  modes  by 
the  cyclic  control  feedbacks. 

2.  The  pilot  acceptance  (rating)  level  is  primarily 

a  function  of  the  automatic  control  and  hover  hold 
functions,  and  does  not  appear  to  be  strongly  de¬ 
pendent  on  the  gross  weight  configuration  of  the 
augmented  stability  aircraft. 

3.  The  best  hover  performance  (position  hold)  capabil¬ 
ities,  within  the  scope  of  the  control  augmentation 
types  considered  here,  are  offered  by  the  longi¬ 
tudinal  cyclic  pitch  (LCP)  control  on  the  stick, 
combined  with  the  conventional  roll  control  to  give 
essentially  a  velocity  control  system  in  the  hori¬ 
zontal  axes.  Pilot  work  load,  and  consequently 
pilot  rating,  is  improved  if  inertial  velocity  feed¬ 
back  is  included  in  the  stability  augmentation  system. 
However,  such  systems  are  inherently  sensitive  to  the 
turbulence,  and  their  total  position  hold  capabilities 
will  degrade  with  the  turbulence  level.  An  improve¬ 
ment  over  the  velocity  control  system,  in  case  of 
precise  hover  hold  tasks,  will  be  a  positive  position 
command  system. 

4.  The  least  desirable  and  the  least  effective  control 
capabilities  are  those  which  employ  pure  attitude 
control,  such  as  are  typical  of  the  differential 
collective  pitch  controlled  tandem-rotor  helicopters 
or  conventional  single-rotor  aircraft. 

5.  The  simulator  model  (Reference  5),  which  is  presently 
used  for  HLH  prototype  work,  has  good  response  fidelity 
at  all  normal  operating  flight  conditions,  as  evidenced 
by  the  response  comparison  with  Model  347  flight  re¬ 
sponse  data.  The  model,  including  the  simplified 
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representation  of  the  external  load  equations 
(Reference  5) ,  should  be  entirely  satisfactory 
for  the  simulator  evaluation  tasks  of  this 
contract. 

6.  The  results  for  the  piloted  simulator  evaluation 
of  the  hover  performance  task  in  general  confirm 
the  validity  of  the  analytical  approach  for  de¬ 
fining  the  relative  performance  capability  of 
various  control  augmentation  modes.  In  absolute 
performance  terms  (RMS,  or  standard  deviation  of 
the  aircraft  position  error) ,  the  analysis  pre¬ 
dicts  well  the  best  pilot  performance  capability 
for  any  specific  augmentation  system.  On  the 
simulator,  the  real  pilot  will  approach  this  per¬ 
formance  level  only  with  sufficient  practice. 

On  the  averaqe,  the  pilot  seems  to  choose  a  per¬ 
formance  point  in  which  the  total  (pilot-aircraft) 
system  damping  level  is  nearly  constant  with  pilot 
gain,  rather  than  the  "optimal"  point  (defined  by 
0.35  damping  ratio  at  minimum  position  error) . 


* 


185 


LITERATURE  CITED 


1.  Seckel,  Edward,  STABILITY  AND  CONTROL  OF  AIRPLANES  AND 
HELICOPTERS,  Academic  Press,  New  York,  1964. 

2.  LONGITUDINAL  HANDLING  QUALITIES  FOR  HOVERING,  Report 

No.  594,  Department  of  Aeronautical  Engineering,  Princeton 
University,  Princeton,  N.J. 

3.  A  NEW  APPROACH  TO  THE  SPECIFICATION  AND  EVALUATION  OF 
FLYING  QUALITIES,  AFFDL-TR-6 9-120 . 

4.  PILOT  DYNAMICS  AND  AIRCRAFT  HANDLING  QUALITIES,  AIAA 
Professional  Study  Series. 

5.  Cogan,  C.,  Gajkowski,  B.  J.,  and  Garnett,  T.  S.,  FULL 
FLIGHT  ENVELOPE  MATH  MODEL  FOR  347/HLH  CONTROL  SYSTEM 
ANALYSIS,  Boeing  Vertol  Report  D301-10148-1,  July  1972. 

6.  A  METHOD  FOR  OBTAINING  AIRCRAFT  MODES  SHAPES  AND  STA¬ 
BILITY  CHARACTERISTICS  WITH  AN  APPLICATION  TO  TANDEM 
ROTOR  HELICOPTERS,  Boeing  Vertol  Report  Aero  inv. 

III-304 ,  January  1966. 


186 


APPENDIX  I 

FULL-FORCE  EXTERNAL  LOAD  MODEL 


ELASTIC  AND  INEXTENSIBLE  CABLE  FORMULATION  FOR 
A  TWO-POINT  SUSPENSION  SYSTEM 


The  method  presented  here  enables  the  evaluation  of  helicopter 
and  external  load  accelerations  at  any  instant  in  time,  given 


the  following  information: 
•  •  • 

XBH'  YBH'  ZBH 

PH'  qH'  fH 


Helicopter  body  axis  velocities 

Helicopter  body  axis  rotational 
rates 


»H'  V  *H 


Helicopter  euler  orientation 
angles 


External  load  body  axis  veloci¬ 
ties 


PL'  qL ' 


External  load  body  axis  rota¬ 
tional  rates 


4» 


L 


f 


External  load  euler  orientation 
angles 


Yl'  Y2 


Forward  cable  angles  relative 
to  inertia  axes 


Y2'  Y4 


Aft  cable  angle  relative  to 
inertia  axes 


Q!s  Forces  and  moments  on  helicopter 

and  external  load,  including 
gravity 


The  following  is  a  list  of  physical  constants  used  in  the  cal¬ 
culations  : 


m 


H 


Mass  helicopter 


ixxH,  iyyh,  izzH,  ixzH 


Helicopter  inertias  relative  to 
body  axes 


m 


L 


Mass  of  external  load 


ixxL,  lYYL/  izzL*  ixzL 


External  load  inertias  relative 
to  body  axes 
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hxl'  hYl' 

hZi 

Location  of  forward  hook  in  heli¬ 
copter  body  axis  coordinates 

hx2'  h*2' 

hZ2 

Location  of  aft  hook  in  helicop¬ 
ter  body  axis  coordinates 

exl'  eYl' 

eZi 

Location  of  forward  sling  and 
load  attachment  point  in  external 
load  body  axis  coordinates 

eX2'  eY2/ 

ez2 

Location  of  aft  sling  and  load 
attachment  point  in  external 
load  body  axis  coordinates 

L1 

Length  of  forward  sling 

l2 

Length  of  aft  sling 

DERIVATION  OF  EXTERNAL  LOAD  EQUATIONS 
Inextensible  Cables 

If  a  dynamical  system  can  be  described  by  a  Lagrangian  L  (qjn, 
qm,  t)  involving  f  coordinates,  and  if  k  constraints  $s  (qi, 
<I2»  •  •  •»  qm»  t)  =  0  are  imposed,  then  the  Lagrangian  equa¬ 
tions  of  motion  can  be  written  as 

* Ufc)  •(%)“*«•  I, >.'*■  (Si'l 

The  Xs(t)'s  are  unknown  functions  of  t,  usually  referred  to  as 
Lagrangian  multipliers.  The  f  equations  of  motion  described 
above,  plus  the  k  constraint  equations,  provide  a  system  of 
(f  +  k)  equations  in  (f  +  k)  unknowns.  The  constraint  equa¬ 
tions  can  be  written  in  dynamical  form  as 


For  the  system  at  hand,  we  chose  the  normal  six-degree-of- 
freedom  coordinates  for  the  helicopter  and  external  load. 
These  plus  two  constraint  equations  representing  the  con¬ 
straints  imposed  by  the  two  slings  result  in  a  system  of  14 
equations . 
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Figure  125.  Full  Force  External  Load  Model. 


Writing  the  normal  six-degree-of-freedom  equations  for  the 
rigid-body  helicopter  as 


mH  XBH  “ 

°XBH  +  all  X1  +  al2  X2 

(1) 

mH  YBH 

°YBH  +  a21  X1  +  a22  X2 

(2) 

mH  ZBH  = 

°ZBH  +  a31  X1  +  a32  X 2 

(3) 

IXXH  Ph 

IxzH  =  °pH  +  a41  X1  +  a42  X2 

(4) 

=  °qH  +  a51  X1  +  a52  X 2 

(5) 

IZZh 

IxzH  P  =  QpH  +  agl  A x  +  ag2  A2 

(6) 

and  a  similar  set  of  equations  for  the  external  load  as 


mL  XBL  = 

QXBl  + 

a71 

X1  + 

a72  X 2 

(7) 

"L  V  - 

qybl  + 

a8l 

X1  + 

a82  X2 

(8) 

mL  5BL  “ 

qzbl  + 

a91 

X1  + 

a92  X  2 

(9) 

IXXL 

-  W* 

’  -  C 

!'  + 
PL 

a101  X1  +  a102  X2 

(10) 

iyyl 

=  Q1  + 

«3L 

alll 

.  X1  +  a112  X2 

(ID 

IZZL 

IxzL  Pl  " 

QK 4 

a121  X1  +  a122  X2 

(12) 

the  physical  significance  of  Ai  and  A?  is  that  they  represent 
the  tensions  in  the  forward  ana  aft  slings  exerted  through  the 
constraint  equations.  The  aj.j's  represent  the  direction 
cosines  of  the  cables  and  relate  directly  to  the  two  con¬ 
straint  equations.  To  develop  these  two  constraint  equations, 
we  observe  the  relative  motion  of  the  helicopter  hook  and  ex¬ 
ternal  load  attachment  point.  Expressing  this  relative  motion 
between  forward  hook  and  forward  attachment  point  in  inertia 
coordinates,  using  rigid-body  motions  of  the  helicopter  and 
load,  we  obtain 
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*xii  "  Axxl  VXL1  +  AXYL  vyLi  +  Axz  VZL1 

"  ^Xh  vXHi  "  axyh  vyHi  ■  axzh  vzHl 

*YIl  =  AYXl  VXLi  +  AYYL  VYLi  +  AYZL  VZLl 

'  ^Xh  vXHi  ■  ayyh  vyHi  "  ayzh  vzHl 


(13) 


(14) 


Rzi!  "  azxl  vxlx  +  azyl  VyLi  +  AzzL  vzLl 


-A  V  -A  V  -A  V 
AZXh  VXHi  “ZYH  VYHi  ZZH  ZHl 


(15) 


We  can  also  describe  this  relative  motion  in  terms  of  two 
angles  and  the  distance  between  hook  and  attachment  point. 

If  we  use  as  the  vector  the  straight  line  joining  the  hook 
and  the  attachment,  with  the  direction  pointing  toward  the 
attachment  point,  then  Y3  will  be  the  angle  through  which  this 
vector  is  rotated  in  order  that  it  lies  in  the  X-Z  inertia 
plan,  and  Yl  will  be  the  angle  to  rotate  this  X-Z  projection 
onto  the  Z  inertia  axis.  Finally,  will  represent  the  mag¬ 
nitude  of  the  vector. 

With  these  definitions,  we  can  write  the  components  (in  inertia 


axes)  of 

this 

vector 

as 

*XIl 

=  L1 

COS  Y3 

sin 

Y1 

(16) 

RYl! 

-  L1 

sin  y3 

(17) 

"m 

=  L1 

COS  Y3 

cos 

Y1 

(18) 

NOTE:  Positive  yi  and  Y3  are  such  that  the  attachment  point 
is  to  the  right  and  forward  of  the  hook. 

We  can  now  differentiate  this  vector  to  yield 

cos  y3  sin  Y^  cos  Y3  cos  Y]_ 

-  Lx  y3  sin  y3  sin  y1  (19) 
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(20) 


RYIi  "  sin  Y3  +  Li  Y3  cos  Y3 
Kl1  =  *1  COS  y3  COS  y1  "  L1  Y1  cos  Y 3  sin  Yl 

-  L1  y3  sin  y3  cos  y1  (21) 

Solving  this  set  of  equations  for  y^,  Y3r  anil  ,  we  get 
I«1  *  Rxix  cos  y3  sin  yx  +  Ryj  sin  y3 

+  RZIi  cos  y3  cos  yx  (22) 

y3  *  ~  **XI;/Ll*sin  y3  sin  Y]_  +  (^yI;l/li)  cos  Y3 

-  (RZJ  /I^)  sin  Y3  cos  y1  (23) 

Y1  "  cos  Y1  “  ^ZIX  sin  Yl)/Ll  cos  y3  (24) 


Substituting  in  our  previous  expressions  for  RXiw  Ryii  and 
Rzij^  in  terms  of  helicopter  and  external  load  motions,  we  now 
have  the  relative  motion  of  the  forward  hook  and  forward  load 
attachment  point  completely  defined.  To  arrive  at  the  first 
constraint  equation,  we  assume  that  =  0  (inelastic  cable)  . 
Using  this  equation,  we  can  arrive  at  our  a.  .'s  associated 
with  the  forward  cable. 

all  "  -^Xh  cosY3sinYl  -  AYXh  sinY3  -  cosy3cosYi 

a21  =  ~AXYh  cosY3sinYl  -  A^  sinY3  -  AZYh  cosy3cosYi 

a31  *  _AXZh  COSY3sinYl  "  ^Zh  SinY3  “  AZZH  COSY3COSYi 

a41  *  _hZl  a21  +  h*l  a3l 

a51  =  hzl  all  _  hxl  a31 
a61  "  "hYl  all  +  hxl  a21 
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71  " 

AXXLCOSY3sinY1 

+  AYXLsinY3 

+  AZXLCOsy3cosy1 

81  " 

AXYLcosY3si”Y1 

+  AY^sinY3 

+  AZYLCOSY3COSYl 

91  " 

AXZLcosY3sinYl 

+  A^sinY-j 

+  AZZiCOsy3cosy1 

101“ 

"eZl  a81  +  eYl 

a91 

111" 

ezl  a71  +  eXl  < 

*91 

121= 

-eYi  a?1  +  eXl 

a81 

We  can  perform  a  similar  analysis  on  the  aft  hook  and  aft 
attachment  point.  The  results  are  similar,  with  the  sub¬ 
scripts  (1)  and  (3)  replaced  by  (2)  and  (4),  respectively. 
Thus,  the  a^j's  associated  with  the  aft  cable  are 

a12  =  “AXXHcosY4siny2  -  AYXHsiny4  -  AZxHcosy4cosy2 

a22  °  "AXYHCOSY4sinY2  "  AYYHsinY4  "  AZYHcosy4cosy2 

a32  =  -AXZHcosy4siny2  -  AYZHsiny4  -  AZZHcosy4cosy2 


42 


— hz«  +  hy5  a. 


l52  “  hz2  a12  “  hxl  a32 


2  “22  ’  x2  “32 
•2  a12 

a62  "  ’hy2  a12  +  hx2  a22 

a72  =  AXXLcosY4siny2  +  AYXLsiny4  +  AZxLcosy4cosy2 
a82  *  AXYLcosy4siny2  +  AYYLsiny4  +  AZYlcosy4cosy2 

a92  "  AXZLCOSY4siny2  +  AYZLsinY4  +  AZZLcosY4COSy2 


(25) 
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a102  ®z2  “82  +  ®Y2  a92 


a112  “  ez2  a72  "  ®x2  a92 

a122  “  "®y2  a72  +  ®x2  a82  {26 

In  the  above  expressions  for  the  aij's,  the  A's  represent 
components  of  euler  transformations  from  body  axes  to  inertia 
axes.  They,  along  with  the  V's,  are  defined  as  follows: 

aXXl  ■  Cos0L  cos\pL 

AXyL  ■  Sin4>L  sin0L  cosi|»L  -  cos<j>L  sin^L 

Axzl  ■  Cos4>L  sin@L  cos^L  +  sin<t»L  sin^L 

AyxL  ■  Cos6l  siniJ>L 

AYyl  =  Sin4>L  sin0L  sinij^  +  cos<f>L  cosij^ 

AyzL  °  Cos4>l  sin9L  sin^L  -  sin#L  cosij>L 

Azxl  =  -sin0L 

AZYl  -  Sin<J>L  cos0L 

AZZL  ■  COS(j)L  COS0L 

AXXH  "  Cos0H  COS^H 

AXYh  =  Sin4>H  sin0H  cos^H  -  cos(j>H  sini|>H 

AXZH  =  Cos*h  sin0H  cos*h  +  sin^H  sin^H 

AyxH  =  C°S0H  sin^H 

AYYh  «  Sin<j>H  sin0H  sin¥H  +  cos<f>H  cosi|»H 

AYZh  *  Cos<j>H  sin@H  sin^H  -  sin<J>H  cosi|>H 
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AzXh  -  -sineH 
AZyH  -  Sin<j>H  cos0H 

AZZh  =  Cos^H  cos0H  (27) 

VXLi  =  XBL  +  eZl  qL  -  eYl  rL 

vYLi  “  YBL  +  exi  rL  "  eZi  PL 

VzLl  “  2BL  +  e*l  PL  -  eXl  qL 

VXL2  =  XfiL  +  eZ2  <3L  “  e*2  rL 

VYL2  "  *BL  +  eX2  rL  -  eZ2  PL 

VZL2  =  ZBL  +  e*2  PL  "  ®X2  qL 

VXH1  *  Xbh  +  hzl  qH  "  h*l  rH 

vyhi  -  ybh  +  hxi  rH  -  hZl  ph 

VzHl  =  ^BH  +  hyl  PH  ~  hXl 

VXH2  =  *BH  +  hz2  qH  ”  h*2  PH 

VYH2  =  YBH  +  hX2  rH  "  hz2  PH 

vZH2  -  ZBH  +  hY2  PH  "  hX2  qH  (28) 

The  two  constraint  equations  are  in  the  form 

**1  “  0  all*BH  +  a21*BH  +  •  •  •  +  a6irH  +  a71xBL 

+  a81*BL  +  *  *  ‘  +  al2irL  <29> 

^2  “  0  “  a12^BH  +  a22^BH  +  •  •  •  +  a62rH  +  a72*BL 

+  a82^BL  +  •  •  •  +  ai22rL  (3°) 
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To  pat  these  constraint  equations  in  the  form  useful  in  solv¬ 
ing  the  system,  we  differentiate  them: 


where 


•  • 

allXBH  +  a21YBH  +  •  *  * 
•  • 

a12XBH  +  a22YBH  +  *  *  * 


•  •  •  • 

QXX  =  "  allXBH  "  a21YBH  *  *  * 
•  •  •  • 

°X2  =  "  a12XBH  "  a22YBH  *  *  * 


Now  we  are  ready  to  solve  the  system  of  equations  for  helicop¬ 
ter  accelerations,  load  accelerations,  and  our  two  unknown 
A's.  First,  we  modify  our  roll  and  yaw  acceleration  equations. 
These  become 

2 

^XXniZZflfiXZjjJPn  =  QpH  +  a4!Xl  +  a42X2 
(IXXHiZZH“IXZH^H  =  °fH  +  a61Xl  +  a62x2 

(IXXLIZZL_IXZL)PL  *  QpL  +  aio,lXl  +  aio , 2X 2 

dxxLIZZL-lxZL) =  °rT  +  a12,lXl  +  a12,2X2 

li 

where 


IzZHQPh  +  IxZHQrH 


+  IxZHQPH 


a41  =  IzZHa41  +  ^H361 
a42  =  IZZHa42  +  IXZHa62 
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a6i  =  IxxHaei  +  IxzHa4i 

a62  =  IxxHa62  +  IXZHa42 

°PL  "  IZZLQPL  +  IxZLQfL 

Qf'L  “  IxxLQrL  +  xxzLQpL 

aio,l  =  IZZLa10,l  +  IXZLa12,l 
a10,2  =  IZZLa10,2  +  rxZLa12 , 2 
a12 , 1  =  IXXLa12,l  +  ^^10,1 
a12,2  =  IXXLa12,2  +  IXZLa10,2 


Our  system  of  equations  is  now  in  a  form  whereby  the  applica¬ 
tion  of  an  elimination-type  algorithm  will  yield  two  simul¬ 
taneous  equations  in  X1  and  \2 • 

all  Xl  +  al2  X2  =  F1 


a21  X1  +  a22  X2  F2 


(32) 


where 


a  =  a 
11  11 


+  a 

H  X1L 


a 


11 


H 


(all  +  a21  +  a31>  .  <a41a41  +  a61a61>  .  a 

-I?.  I 


m 


H 


WZZH  ~XZH 


51 

yyh 
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I 


a 


Hi 


(a71+a81+a91)  .  {a10,la10,l+a12,la12,l)  ,  all,l 
-  +  - i - 1 - -  +  i - 

XXL  ZZL  XZL  YYL 


"L 


a12  “  a12 


H 


+  ot 


12, 


a 


12 


(alla12+a21a22+a31a32)  .  (a41a42+a61a62)  ,  a51a52 

'  '  A  *  ~  ’ 


H 


m 


H 


I  I  -i2 
XXH  ZZH  XZH 


YY 


H 


(a71a72+a81a82+a91a92)  'a10 , la10 , 2+a12 , la12 , 2 [ 

12l  mL  IXXLIZZL"IXZL 


all,lall,2 

jyyt 


F1  =  Qx  "  Fi  ■  Fi 

1  A1  1H  iL 


(ailQXBH+a2lQYBH+a31QZBH1 


(ailiQ~  +aciQr  )  aciQ' 


41  pR  61  r 


H 


51wq 


H 


H 


m 


H 


I  I  -1^ 
xxxHxzzH  XXZH 


YY 


H 


(a71QXBL+a81QYnT+a91CZnT  }  (a10,lQpT+a12,lQrT  } 


BL 


BL 


mT 


I  I  -i^ 
XXL  ZZL  XZL 


an  iQ' 
+  ^ ^ 

IYY, 


a21  =  a21H  +  a21L 


a 


21 


H 


(a12all+a22a21+a32a31)  .  (a42a41+a62a61}  .  a52a51 

-  +  - ^ -  +  — - 


m. 


H 


IxxhIzzh'Ixzh 


YY 


H 
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a 


(a72a71+a82a81+a92a91) 


21. 


m. 


^a10>2a10»l+a12/2a12#l^ 

I  I  -l2 
XXL  ZZL  XZL 


a22  "  a22H  +  a22L 


.  all «  2all , 1 

I 

YYL 


a 


22 


H 


(a12+a22+a32)  .  (a42a42+a62a62)  .  52 

-  +  - - - - — 2 -  +  - 

mH  IXXHIZZH_IXZH  IYYh 


*  A  a  a^i 


,) 


a 


22, 


72  b2  92  x“10,2“10,2  “12,2“12,2'  +  “11,2 

I„„  I. 


m. 


bXX  ~ZZ_  *XZT 

li  Li  Li 


yyt 


F2  =  QX  "  FZ  "  FZ 
2  A2  ZH  ZL 


■H 


(a12QXBH+a22QYBH-a32Q2BH>  <a42QpH+-62Q?H5 
-  + - 


a52°q 


H 


m 


H 


IxXBIzZH  'XZH 


YY 


H 


(a72QxBT+a82QyBT+a92QzBT  *  (a10 , 2QpT **12 1 2QrT * 

JLl  li  Li  .  Jj  Lj 


m. 


I  I  -i^ 
XXL  ZZL  XZL 


aH,2QqL 

IYY. 


The  solution  is 


L1  "  (Fla22  "  a12F2)/(a22all  "  a12a21) 


(33) 


X2  "  (allF2  "  Fla21)/(alla22  “  a12a21) 


(34) 
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Substitution  of  these  values  into  the  appropriate  equations 
will  enable  us  to  solve  for  the  unknown  accelerations. 


Ybh 

*bh 


Ph 

qH 


XBL 


yBl 

*bl 


Pl 

qL 


<QxBh 

+  a11X1  ■ 

f  a12X2^mH 

(QyBh 

+  a21X1  ■ 

*■  a22X2)/mH 

<Qzbh 

+  a31X1  ■ 

f  a32X2^mH 

(q;h  + 

a41Xl  + 

a42X2^  ^IXXhIZZjj  ”  JXZh^ 

(q*h  + 

a51Xl  + 

a52X2)^IYYH 

(QM  + 
H 

a61Xl  + 

a62X2^  ^  ^XXji^ZZjj  "  IXZh^ 

<QxBl 

+  a71X1  +  a72X2)/mL 

o 

K 

a 

t* 

+  a^Xj^  H 

h  a82X2)/mL 

(qZbl 

+  a91Xl  +  a92X2^mL 

<Q»1  + 

a10,lXl 

+  a10,2X2)/{IXXLIZZL  _  ^Zl) 

all,lXl 

+  all,2X2)/;CYYL 

(op  + 

L 

a12,lXl 

+  a12,2X2)/(IXX  IZZ  “  IXZ  J 

Li  Lt  Li 

(35) 


Now,  since  the  helicopter  body  axes  and  also  the  external  load 
body  axes  are  rotating  axes,  the  inputs  into  the  integrators 
to  yield  translational  velocities  are 

*BH  “cT^Bh  "  ^BH  qH  +  YBh  rH*  dt 

Yfijj  =J*(Ybh  ”  rH  +  ZgH  PB)  dt 

Zbh  =  J* ^Bh  ”  YBh  Ph  +  XgH  qH)  dt 
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(36) 


XfiL  =  /(X Bl  -  Zbl  qL  +  yBl  rL)  dt 

*BL  =  /(YbL  '  XfiL  rL  +  ZBL  Pl)  dt 

*fiL  =/  (Z’BL  "  YBL  pl  +  Xbl  qL)  dt 


eralized  forces^Q^sK ^These^hall^I^f ' iTly,d<rfine  the  gen- 
ways:  ese  sha11  be  defined  in  the  following 


Qxb„ 

“  XH 

-  mHg  sin 

6H 

=  yh 

+  mHg  cos 

9H 

sin  4>h 

Qzr 

bh 

=  ZH 

+  m^g  cos 

9H 

cos  # 

XI 

qph 

*  lh 

Q' 

PH 

=  qph 

+  JXZH  PH 

% 

(Izz 

Q9H 

-  mh 

Q<jH  \  "  IxzH  (PH2-rH2)  -  <JXX  -hz  >  rH  PH 

**  H 


‘H 

QrH  *  N« 


QfH  °rH  IxZH  qH  ^  ~  (IyyH~IxxH)  PH  qH 

°XBl  =  XL  “  mL  g  sin  ®L 

QyBl  "  YL  +  ®L  g  COS  eL  sin  $L 

zl  +  mL  9  cos  0^  cos  $L 


°PL  *  Ll 

Pl  °Pl  +  JXZL  PL  qL  “  ^ZZt^YYt*  qL  PL 


Pl  *■*» 

Q9L  *  “l 
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+ 


°qL  ”  Q(2L  +  IxZl  (?L  ”rL  *  ~  (Ixxl"IzZl)  rL  Pl 

Qt,  “  N, 


Qf, 


°rL  "  IxzL  qL  rL  "  (iyyl“ixxl)  pl  q; 


(37) 


where  Xjj,  Yh#  Z y,  Ljj,  Mh,  and  Njj  represent  the  forces  and 
moments  exerted  on  the  helicopter  from  such  things  as  fuselage 
aerodynamics,  rotor  thrust,  and  rotor  hub  moments.  Xl,  Yl, 

ZL,  lL'  Ml,  and  Nl  represent  the  forces  and  moments  exerted 
on  the  external  load  from  such  things  as  load  aerodynamics. 


The  preceding  formulation  is  valid  for  all  motions  of  the 
nelicopter  and  load  if  we  assume  rigid  cables.  This  is  true 
since  the  X 1 s  will  always  be  calculated  to  maintain  the  dis¬ 
tance  from  hook  to  attachment  point  constant.  In  other  words, 
if  there  arose  a  situation  in  which  the  forces  on  the  load 
were  such  that  the  tendency  would  be  for  the  attachment  point 
to  move  closer  to  the  hook,  the  cables,  through  the  X's,  would 
exert  a  compressive  force  to  maintain  the  distance  at  cable 
length.  To  allow  for  flexible  cables,  we  must  eliminate 
constraint  equations  when  cable  forces  become  compressive. 

This  is  done  by  simply  setting  the  X's  to  the  maximum  com¬ 
pressive  force  assumed  for  the  cables  (for  noncompress ive 
cables,  the  X's  are  set  to  zero).  Thus,  any  degree  of  rigid¬ 
ity  of  the  cables  may  be  assumed. 


Thus,  the  problem  of  allowing  the  cables  to  go  slack  is  easily 
handled.  The  problem  of  returning  the  cables  to  tension  when 
the  conditions  call  for  it  is  somewhat  more  complicated.  The 
complication  arises  because  we  cannot  suddenly  begin  the  X's 
calculation  at  some  instant  in  time  without  assuring  that  the 
velocities  of  the  helicopter  and  load  conform  to  the  con¬ 
straint  equations.  To  solve  this  problem,  we  must  reinitialize 
the  helicopter  and  load  velocities  and  rotations  at  the  instant 
the  cables  become  taut.  We  do  this  by  applying  impulse- 
momentum  theory.  The  mathematical  expression  that  we  will 
use  to  relate  the  velocities  before  tautness  to  after  tension 
is  returned  is 


(38) 
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In  this  expression,  (3L/3qj)f  represents  the  momentum  af£er 
impulse  in  the  direction  of  the  qj  variable,  while  (3L/3qj)0 
represents  the  momentum  before  impulse  in  the  qj  direction, 
and  Ij  represents  the  component  of  the  impulse  in  the  qj 
direction. 

In  terms  of  our  specific  problem,  we  have 


(1) 

VCsHf  - 

mHXBHD  “  allIl  +  a12Z2 

(2) 

Vbh,  - 

mH^BHQ  “  a21Il  +  a22I2 

(3) 

mHZBHf  ' 

mH^BHo  "  a31Il  +  a32I2 

(4) 

(I  P 

XXH  Hf 

-i  r  )-(i  p  -i  r  ) 

XZHV  XXH  Ho  XZH  V 

"  a41Zl  + 

(5) 

IyVH  X 

IYYU  qH  =  a51Il  +  a52I2 

(6) 

-i  p  )-(i  r  -i  p  ) 
XVV  Z3H  Ho  XZH  Ho 

=  a61Zl  + 

(7) 

mLXBL£  - 

mLX0Lo  *  a71Zl  +  a72Z2 

(8) 

mLyBLf  * 

mL^BLo  "  a81Il  +  a82I2 

(9) 

",lzbL£  - 

mLZBLQ  =  a91Il  +  a92Z2 

(10) 

_IXZ  rL  PL  _IXZ  rL  * 

L  f  L  o  L  o 

=  a10,lIl 

+  a10,2Z2 

(11) 

jyyl  qf 

ZYYt  qL  "  all,lIl  +  all, 

Ju  O 

2I2 

(1?) 

UzZLrLf 

_IXZ  PLf)-(IZZ  rL  _IXZ  PL  } 
L  f  L  O  L  O 

=  a12 , l1! 

+  a12,2I2  (39) 
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Since  we  know  the  velocities  and  rotations  before  the  impulse, 
*BHo'  *BHo*  •  •  •/  Tl0/  and  if  we  assume  that  the  direction 

coefficients,  an,  a^2»  •  •  •»  ai2  2'  do  not  change  during  the 

impulse,  we  see  that  we  have  12  equations  in  14  unknowns. 
However,  the  final  velocities  and  rotations,  &bh£ •  $BHf» 

.  .  . ,  must  conform  to  the  following  constraint  equations 

(13)  a^Xe^  +  a21^BHf  +  *  *  *  +  all,lqLf  +  a12,irLf  =  0 

U4)  +  a22*BHf  +  •  •  •  +  all,2qLf  +  a12,2rLf  =  0 

(39) 

Thus,  we  have  a  system  of  14  equations  in  14  unknowns  which 
can  readily  be  solved. 

When  the  external  load  is  restrained  within  the  hull  of  a  sta¬ 
tionary  ship,  the  following  equations  are  superimposed: 


\  "  \  “  PL  "  qL  "  rL  "  0 

a71  “  a81  =  a10,l  =  all,l  =  a12 , 1  = 

a72  “  a82  "  a10 , 2  =  all,2  =  a12 , 2  = 


(40) 


Elastic  Cable  Formulation 


The  preceding  derivation  for  inelastic  cables  can  easily  be 
modified  to  handle  elastic  cables.  In  fact,  the  calculations 
are  somewhat  shortened.  The  difference  arises  in  that  for  the 
inelastic  cables,  we  forced  £<l  =  1*2  =  0;  and  by  doing  this, 
we  were  able  to  define  two  additional  system  equations  re¬ 
lating  helicopter  and  load  motions.  With  these  two  additional 
L  equations,  we  were  able  to  calculate  the  necessary  restrain¬ 
ing  tensions  (X^,  Xj)* 

In  the  case  of  elastic  cables,  we  allow  the  distance  between 
hook  and  attachment  point  to  vary,  as  if  helicopter  and  load 
were  two  completely  independent  free  bodies.  We  then  measure 
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the  distance  between  hook  and  attachment  point  and  the  rate  of 
change  of  this  difference,  exactly  as  we  did  for  the  inelastic 
cable  case;  that  is, 


Lx  =  RXll  cos  y3  sin  y1  +  Ry^  sin  y3  +  Rz^  cos  Y3  cos  yx 


L2  =  Rxj  cos  y4  sin  y2  +  RYj;  sin  y4  +  Rzj  cos  Y4  cos  Y2 
2  2  2 


However,  unlike  the  inelastic  case  where  the  X's  were  calcu¬ 
lated  as  necessary  internal  restraining  forces,  the  X's  for 
the  elastic  cable  case  are  calculated  directly  from  L]_,  L2, 
L,  ,  and  £2  an£*  the  physical  properties  of  the  cables.  If  we 
assume  that  the  cables  will  exert  a  force  proportional  to 
change  in  length  and  rate  of  change  in  length,  then  we  can 
calculate  the  X's  in  tension  as 


(1)  X  x  =  -kT  ALX 


(2)  X2  =  -kT  AL2 
2 

and  in  compression  as 


(3)  \1  =  -kc  ALX 


(4)  X2  «  -k^  4L, 


(41) 


where  k^*  ^T2»  and  are  cable  spring  constants; 

Ct2 •  Ccj.,  an<*  Cc2  are  ca^^e  damping  rates. 

Once  these  X's  are  calculated,  the  procedure  is  the  same  as 
for  the  inelastic  case. 

Digital  Computer  Mechanization  of  External  Load  Equations 

The  sling  load  equations  were  programmed  on  the  IBM  360-65. 
Numerical  integration  of  these  helicopter  and  load  accelera¬ 
tions  was  performed  using  fourth-order  Runge-Kutta  methods. 
The  following  is  a  brief  description  of  the  segments  of  this 
program,  indicating  one  possible  sequence  of  calculations. 
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1. 


Main  Routine  -  Executive  program  used  to  read  in  physical 
parameters  of  helicopter  and  load,  set  up  initial  condi¬ 
tions,  set  up  integration  cycle,  and  plot  output. 

2.  Trim  Subroutine  -  Subroutine  called  from  main  program 
during  setting  up  of  initial  conditions. 

3.  TENSNS  Subroutine  -  Subroutine  called  from  main  program 
during  integration  cycle  to  calculate  cable  tensions. 

4.  Impact  Subroutine  -  Subroutine  called  from  main  program 
during  integration  cycle  if  either  or  both  cables  are 
going  from  slack  condition  to  taut  condition.  This  sub¬ 
routine  calculates  helicopter  and  load  velocities  and 
rotations  after  tautness  occurs,  based  on  conditions  just 
before  tautness. 

5.  RKS  Subroutine  -  Subroutine  employing  fourth-order  Runge- 
Kutta  method  to  numerically  integrate  our  equations  of 
motion.  Subroutine  is  called  from  main  program  during 
integration  cycle. 

6.  Plot  Subroutine  -  Subroutine  called  by  main  program  to 
plot  output  data  on-line  on  the  IBM  360-65. 
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APPENDIX  II 

PARAMETRIC  HELICOPTER  DATA 


The  turbulence  is  represented  as  random  perturbations  in  lon¬ 
gitudinal  and  lateral  velocity.  The  response  of  a  linear 
system  to  a  random  input,  in  terms  of  power  spectra,  is 

PQ  (w)  =  |G  (jw)  |  2  Pj  (w) 


where  G{s)  is  the  transfer  function  of  the  system 

P^(w),  P  (w)  are  the  power  spectrums  of  the  input 
and  output,  respectively 


The  mean-square  response  of  the  output  is 

a  ac 

''o'  I?  /  Po  <“>  d“  =  SF  /|G(jo)|2  Pjlo)  do 

—  oc  — oc 


As  an  example,  suppose  it  is  desired  to  calculate  the  mean- 
square  longitudinal  position  error  response  to  a  one¬ 
dimensional  gust,  assuming  the  following  gust  spectrum: 


Pz  (u>) 


2og2  wN 

2  .  2 

U)  + 


where  uN  is  gust  break  frequency 

og  is  the  standard  deviation  of  gust 


The  desired  system  transfer  function  is 

X 

G  (s)  =  =£■  (s) 

ug 
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and  the  resultant  mean-square  position  error  response  is 


_L  f 

2tt  * 


uf 


-  2 
2  ag  wN 

“T~  2 

a)  +  a)N 


da3 


To  perform  this  analysis  and  also  to  determine  the  stability 
of  the  system  in  terms  of  eigenvalues,  phase  margin,  etc.,  a 
digital  program  was  written  with  the  following  options. 

Calculations  of 

(1)  Eigenvalues  and  eigenvectors  of  state  matrix. 

(2)  Zeroes  and  poles  of  desired  transfer  function, 
including  elimination  of  like  poles  and  zeroes. 

(3)  Frequency  response  of  transfer  function. 

(4)  Inverse  Laplace  coefficients  of  transfer  function. 

(5)  Time  response  of  transfer  function  to  random 
disturbance  (turbulence) . 
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APPENDIX  III 

SUMMARY  OF  SIMULATOR  EVALUATION  DATA 


This  appendix  presents  the  simulator  evaluation  run  log  (Table 
XIX)  ,  and  includes  the  list  of  all  runs  which  were  directly 
used  to  obtain  the  quantitative  and  the  qualitative  data  sum¬ 
marized  in  the  main  sections  of  this  report. 
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Date  6-27-73 _  Constants:  A/C  Gr  Wt  50,000  Lb  Notes:  Turbulence  Level  is  Ug  and 

Configuration  50,000  Lb^  External  Load  None  Wg  »  0  For  All  Runs 

Wind  10  KT 
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lable 
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TABLE  XIX.  ATH  SIMULATION  LOG  (Continued) 
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TABLE  XIX.  ATH  SIMULATION  LOG  (Continued) 
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Control 


TABLE  XIX.  ATH  SIMULATION  LOG  (Continued) 


